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”  A '^Electronic  breadbanding  is  a  technique  whereby  very  low  frequency 
(VLF)  antennas  c..u  be  resonated  (rctuned)  synchronously  with  a  f requeue 
shif t-kovrd  (FS::.)  signal  for  transmission  of  binary  data  at  rates  well 
over  1 000  baud.  The  purpose  of  this  thesis  is  to  investigate  the  crans 
minted  energy  density  spectrum  for  an  electronically  broadbanded  VL1 
antenna  system  in  which  a  frequency  shift  can  occur  only  at  the  instant 
of  an  antenna -current  zero-crossing  with  positive  slope,  A  computer 
solution  shows  that  the  transmitted  spectrum  is  approximately  equiva¬ 
lent  to  the  spectrum  of  an  ideal,  constant  amplitude,  FbK.  signal, 

A  further  investigation,  again  by  computer  program,  snows  that,  j.t 
the  total  frequency  shift  between  the  marking  and  spacing  frequencies 
occurs  over  a  finite  period,  rather  than  instantaneously,  cben  the 
spectrum's  sidelobcs  will  be  reduced  significantly.  ihe  spectrum  for 
a  27,000  Hz  center-frequency  VLF  antenna  transmitting  at  1600  baud  by 
means  of  electronic  broadbanding  can  be  reduced  from  3  to  18  dB  in  the 
first  five  sidelobcs.  No  evaluation  is  made  on  the  detectability  of 
the  signal  with  the  non- instantaneous  frequency  shift.  Numerous  plots 
of  the  computer-generated  spectra  are  included, 
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Electronic  broadband Ing  is  a  technique  whereby  very  low 
frequency  (VLF)  antennas  can  be  resonated  (retunod)  synchro¬ 
nously  with  a  frequency -shift -Keyed  (FSk)  signal  for  trans¬ 
mission  of  binary  data  at  rates  well  over  1000  baud.  The 
purposo  of  this  thesis  is  to  investigate  the  transmitted 
energy  density  spectrum  tor  an  electronically  broadbanded 
VLF  antenna  system  in  which  a  frequency  shift  can  occur  only 
at  the  instant  of  an  antenna -current  zero-crossing  with  pos¬ 
itive  slope.  A  computer  solution  shows  that  the  transmitted 
spectrum  is  approximately  equivalent  to  the  spectrum  of  nn 
ideal,  constant  amplitude,  FSK  signal. 

A  further  investigation,  again  by  computer  program, 
shows  that  if  the  total  frequency  shift  between  the  marking 
and  spacing  frequencies  occurs  over  a  finite  period,  rather 
than  instantaneously,  then  the  spectrum's  stdelobes  will  be 
reduced  significantly.  The  spectrum  for  a  27,000  Hz 
center -frequency  VLF  antenna  transmitting  at  1600  baud  by 
means  of  electronic  broadbanding  can  be  reduced  from  3  to  IS 
dB  in  the  first  five  siuelobes.  No  evaluation  is  made  on 
the  detectability  of  the  signal  with  the  non-instantanecua 
frequency  shift.  Numerous  plots  of  the  computer-generated 
spectra  are  included. 
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I ,  I nt  roductl  on 


The  United  States  Air  Force's  Survivablc  Low  Frequency 
Communication  System  uses  the  very  low  frequency  (VLF)  por¬ 
tion  of  th"  spectrum  to  achieve  reliable  lone-ranee  commu¬ 
nications.  Ihis  thesis  is  an  examination  of  the  energy 
spectrum  of  a  specific  FSR  modulation  scheme  that  could  be 
used  for  VLF  communications.  In  this  thesis,  frequency- 
shift  Keying  is  defined  as  a  form  of  binary  frequency  modu¬ 
lation  in  which  a  specific  frequency  is  transmitted  for  a 
digital  zero,  and  a  different  frequency  is  transmitted  for 
a  digital  one. 

In  1972,  Home  Air  Development  Center  sponsored  a  devel¬ 
opment  contract  for  the  purpose  of  building  a  prototype 
mechanism  that  could  be  used  to  broaden  the  apparent  band¬ 
width  of  conventional  VLF  antennas  so  that  frequency- shift- 
keyed  CFSK)  signals  could  be  transmitted  at  higher  rates 
than  those  previously  possible.  The  mechanism  was  :  ?..t.c2  ui« 
Experimental  Transmitting  Antenna  Modulator  (ETAM)  and  its 
design  was  made  possible  by  the  use  of  newly-developed , 
high-speed,  silicon-controlled  rectifiers  (SCR).  The  ETAM 
Is  a  realization  of  a  syncinonous-reeonation  broad banding 
technique  wherein  the  antenna  is  resonated  (retuned)  in 
accordance  with  the  frequency  o£  the  FS'U  signal  being  trans¬ 
mitted.  In  addition,  each  frequency  change  can  occur  only 
at  the  instant  ol  an  antenna-current  zero-crossing  with 
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positive  slope  (Ref  4 i 1-15),  The  purpose  of  this  thesis  is 
to  determine  the  frequency  spectrum  for  FSK  signals  trans¬ 
mitted  from  a  VLF  antenna  system  which  uses  an  electronic 
broad banding  device  such  as  the  EXAM. 

The  study  is  limited  to  VLF  transmitting  systems  in 
which  there  is  no  amplitude  modulation  of  the  transmitter 
voltage  envelope.  However,  in  contrast  to  FSK  signals  whose 
entire  frequency  shift  occurs  instantaneously,  this  study 
Includes  an  analysis  of  electronically  broadbanded  systems 
whose  frequency  transitions  are  non- Instantaneous, 


II ,  Background 


One  method  of  modulating  radio  signals  for  transmission 
of  digital  data  is  frequency -shift  keying,  wherein  a  spe¬ 
cific  frequency  is  transmitted  for  a  mark  (a  digital  one) 
and  a  different  frequency  is  transmitted  to  signify  a  space 
(a  digital  zero).  Transmission  rates  for  FSK  signals  in  the 
VLF  range  (10  -  100  kHz)  are  severely  limited  by  the  band¬ 
width  of  conventional  VLF  antennas.  Because  of  the  extreme¬ 
ly  long  wave  lengths  (up  to  3C  km),  VLF  antennas  are  electri¬ 
cally  very  short,  and  in  order  to  achieve  efficient  power 
transmission,  the  antenna/transmitter  system  must  be  tuned 
to  form  a  high-Q  resonant  circuit.  The  result  is  that 
typical  VLF  antennas  may  have  a  3  dB  bandwidth  of  100  Hz  or 
less  (Ref  b$33).  The  bandwidth  limitation  restricts  the  FSK. 
transmission  rate  to  approximately  50  baud,  or  to  a  total  of 
50  marks  and  spaces  per  second,  because  of  the  phase  distor¬ 
tion  and  amplitude  transients  which  would  occur  at  higher 
transmission  rates  (Ref  4i555). 

Numerous  methods  have  been  investigated  to  increase 
VLF  transmitting  system  bandwidths,  but  the  most  successful 
technique  is  to  retune  the  antenna  system  in  synchroni¬ 
zation  with  the  FSK  signal  by  either  capacitance  modulation 
or  by  inductance  modulation  (Ref  4*555),  The  method  of  syn¬ 
chronous  retuning,  or  synchronous  resonation,  can  be 
described  in  terms  of  a  mathematical  model. 
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A  VLF  Antenna  Model 

VLF  antennas  can  be  modeled  as  series  RLC  circuits, 
and  a  generally  accepted  model  i3  shown  in  Fig.  1.  The 
model  is  valid  for  frequencies  well  below  the  uncompensated 
resonance  frequency  of  the  antenna.  Due  to  the  effects  of 
the  electromagnetic  fields,  the  antenna's  apparent  in¬ 
ductance  varies  with  frequency.  For  frequencies  which  are 
within  the  VLF  range  and  that  are  less  than  half  the  uncom¬ 
pensated  resonance  frequency  of  the  antenna  (approximately 
100  kHz),  the  apparent  inductance  changes  by  less  than  5% 
and  is  considered  as  a  constant  (Ref  7 i 19,645), 


R. 


R. 


R. 


R. 


R 


g 


C 


transmitter  internal  resistance 
copper-loss  resistance 
radiation  resistance 
ground-loss  resistance 
apparent  antenna  inductance 
tuning  inductance 
antenna  static  capacitance 


Fig.  1,  Very-Low-Frequency  Antenna  Model 
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The  average  total  power  (Pt)  transmitted  from  a  VLF 
antenna  can  be  expressed  as 

Pt  K  40  tt2  I2  (h/X)2  watts  (1) 

where  I  is  the  antenna  current  peak  value,  h  is  the  effec¬ 
tive  height  of  the  antenna,  and  X  is  the  transmitted  wave¬ 
length,  From  Eq  (1),  a  radiation  resistance  (R  )  may  be 
defined  as  the  resistance  which  would  dissipate  P  watts 
when  a  sinusoidal  current  with  peak  value  I  flows  through 
it  (Ref  7i645).  Hence, 

Rr  -  (2Pt)/l^  ■  80  it2  (h/X)2  ohms  (2) 

The  radiation  resistance  for  installed  VLF  antennas 
is  normally  between  0,05  and  0.3  ohms,  while  the  total  of 
all  the  other  circuit  resistances  (Rrnl_)  is  between  0.1 
and  0.7  ohms.  The  result  is  that  the  transmission  effi¬ 
ciency  (Rr/p-tot)  is  in  the  0,25-to-0.65  range.  However, 
one  modern  VLF  antenna  has  a  transmission  efficiency  of 
0.88  (Ref  8 « 142). 

The  tuning  inductance,  Lt,  is  added  to  the  circuit 
in  order  to  lower  the  system  resonant  frequency  from  its 
uncompensated  one  to  the  desired  frequency  in  the  VLF  range. 
Such  inductors  carry  large  currents  and  have  been  conven¬ 
tionally  constructed  as  air -core  coils. 

Since  the  circuit  is  a  familiar  RLC  circuit,  the 
system  3  dB  Bandwidth  is  the  ratio  of  the  resonant  fre¬ 
quency  to  the  system  Q.  The  system  Q  is  the  ratio  of  the 
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circuit  reactance  to  the  total  resistance.  The  equation  for 
the  bandwidth  is 

Bandw*idth(3dB)  «=  2tr  f2  C  RCot  Hz  (3) 

References  to  bandwidth  are  to  be  made  with  caution, 
since  it  must  be  explicitly  stated  whether  the  bandwidth 
referred  to  is  the  system  bandwidth,  or  that  o1  the 
antenna  alone.  The  system  bandwidth  is  twice  the  antenna 
bandwidth,  due  to  the  transmitter  being  matched  to  the 
antenna,  which  appears  as  a  resistance  at  the  system 
resonant  frequency.  In  this  paper,  the  system  bandwidth 
will  be  used  unless  otherwise  stated, 

The  narrow  bandwidth  implies  that  for  FSK  modulation, 
the  mark  and  space  frequencies  must  both  be  well  inside 
the  3  dB  bandwidth  limits,  and  that  the  baud  rate  must  be 
low  enough  so  that  the  transmitted  signal  amplitude  and 
phase  are  not  unduly  distorted  by  transients.  Such  tran¬ 
sient  distortion  reduces  radiated  power  and  receiver 
detection  efficiency  (Ref  /ti556).  One  method  for  increasing 
the  apparent  bandwidth,  and  thereby  increasing  the  potential 
baud  rave,  is  that  of  synchronous  resonat-ion. 

Energy  Relationships  for  Synchronous  Resonation 

Higher  transmission  rates  and  wider  apparent  bandwidths 
may  be  achieved  for  pulsed  transmission  circuits,  such  as 
for  FSK  modulation,  by  a  variation  of  the  following  concept. 
When  it  is  possible  to  insert  or  to  remove  a  portion  of  one 
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type  of  reactive  element  with  the  appropriate  amount  of 
stored  energy  at  instants  when  the  other  type  of  reactive 
element  in  a  resonant  circuit  has  dissipated  all  of  its 
stored  energy,  then  an  instantaneous  steady- state  oscilla¬ 
tion  is  achievable.  Such  a  technique  permits  instantaneous 
switching  between  a  quiescent  condition  and  a  steady- state 
response.  Further,  it  enables  a  change  in  resonant  fre¬ 
quency  of  an  oscillating  circuit  as  often  as  once  each  cycle 
In  the  application  to  FSK.  transmission,  the  transmitter  fre¬ 
quency  and  the  circuit's  resonant  frequency  may  be  simul¬ 
taneously  changed  without  exciting  large  transients,  by 
inserting  or  removing  a  portion  of  the  reactive  element  whoa 
stored  energy  is  at  a  minimum.  The  result  is  that  the 
steady- state  response  is  obtained  almost  immediately  (Ref  2i 
375),  The  physics  of  such  a  switching  operation  is  apparent 
when  viewed  in  a  set  of  differential  equations, 

A  circuit  model  from  which  the  circuit  differential 
equation  can  be  derived  is  shown  in  Fig.  2.  Each  circuit 
has  a  switchable  DC  source  (V  or  I)  for  establishing  initial 
conditions,  but  both  circuits  have  the  same  differential 
equation  when  their  respective  switches  are  in  position  2. 
The  equation  is 


v(t) 


L  ^ 
L  dt 


+  i  R  + 


hi: 


(1)  dt 


(4a) 


or. 


v(t) 


BLq+Rq4  (q. 


+  q) 


(4b) 


where  1  is  the  time-dependent  current  and  q  is  the  time- 
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VJ* 


dependent  capacitor  charge.  If  t  or  q  i s  assumed  to 
exhibit  a  zero -transient  response,  then  the  required 
initial  conditions  can  be  determined  from  Eqs  (4a)  and  (4b). 

Assume  for  Case  A,  the  switched-C  circuit  of  Fig.  2, 
that  the  steady-state  current  is  ia(t)  «=  losin(wt)  and  that 
the  switch  SWl  is  moved  from  position  l  to  position  2  at 
time  t  =  0,  For  Case  8,  the  switched-L  circuit  of  Fig.  2, 
assume  the  steady-state  current  is  ib(t)  =  Iocos(wt),  and 
that  the  switching  is  the  same  as  for  the  switched-C  circuit. 

The  result  for  Case  A  is 

v_(t)  *  wLI^  cos(wt)  +  RI_  sin(wt) 

O  U 

+  £§  (l-cos(wt))  +  ^  J  ia(t)  dt  (5) 
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and  the  result  for  Case  B  is 

vb(t)  *=  -wLIQ  sin(wt)  +  RIq  cos(wt) 

+  ^  sin(wt)  t  j*  f  ib(t)  dt  (b) 

-OP 

For  a  resonant  cLrcuit,  wL  =  which  enables  cancellation 
of  two  terras  in  each  of  Eqs  (5)  and  (6).  The  results  are 

va(t)  ■=  RIq  sin(wt)  +  ^0  +  £  f  ia(t)  dt  (7) 

xJ  «.  CO 

and 

i  r° 

vb(t)  «=  R1  cos(wt)  -f  g  I  i5(t)  dt  (8) 

J  -00 

for  Cases  A  and  B,  respectively.  The  integral  terms  of 
Eqs  (7)  and  (8)  are  the  initial  voltages  on  the  capacitors. 

If  the  voltage  generator  V  for  Case  A  equaled 

v  -  h  fVc>  dt  -  -  k  <9> 

+J  »oo 

and  if  the  initial  current  i_(0)  **  0,  or  if  the  capacitor 

ft 

voltage  for  Case  B  was 

i  f  °i5(t)  dt  ■  0  (10) 

and  the  current  source  I  «  ib(0)  E  IQ.  then  ia(t)  and  ib(t) 

would  be  zero-transient  sinusoids.  If  different  assumptions 

had  been  made,  such  as  ia(t)  **  IQcos(wt)  or  i(b)  «*  lQsin(wt), 

then  the  exact  initial  conditions  I  =  i_(0)  *=  1  (for  the 

i„  a  0 

switched-C  case)  and  V  *=  -  ^  (for  the  switched-L  case) 
required  for  a  zero -transient  response  could  not  be  obtained, 
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The  application  of  the  above  principle  to  a  circuit 
operating  at  resonance  can  result  in  a  frequency-change 
current  transient  that  is  a  small  fraction  of  the  steady- 
state  current  amplitude.  Fig.  3  depicts  the  relations 
which  yield  such  a  shift.  To  make  a  transient-less  fre¬ 
quency  change  to  IQsin(wt) ,  the  current  must  be  zero  and 
the  capacitor  voltage  must  be  -  Such  constraints  on  the 

current  are  satisfied  every  2ivrr  radians,  but  the  capacitor 


Fig.  3.  Steady-State  Current  Relations 

voltage  constraint  is  not  met  exactly.  The  capacitor  is 
charged  to  -I0/(w0C),  where  wQ  is  the  radian  frequency  prior 
to  the  change,  and  such  points  are  designated  by  "A”  on  Fig. 
3,  To  make  a  change  to  IQcos(wt),  the  capacitor  voltage 
must  be  zero  and  the  current  must  be  at  a  maximum.  Such 
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time  points  are  designated  as  "B”  on  Fig.  3,  A  change 
to  I0cos(wt)  would  have  no  transient  because  the  current 
peak  is  independent  of  the  resonant  frequency.  Therefore, 
an  ideal  frequency  change  can  be  made  by  modulating  the 
capacitor  when  its  voltage  is  zero,  or  a  near-ideal  fre¬ 
quency  change  can  be  made  by  adjusting  the  inductor  when 
the  current  is  zero  and  the  capacitor  V0iL<i6e  is  at  ~ 

negative  peak  (a  minimum).  For  the  two  cases  discussed 

2 

in  this  paragraph,  the  capacitor  energy  (CV  /2)  and  the 

2 

inductor  energy  (LI  /2)  are  at  mutually  exclusive  maxima, 
that  is,  one  or  the  other  has  all  the  stored  energy  at  the 
time  of  the  change  (Ref  2i376). 

Synchronous  resonation  schemes  have  been  attempted 
both  through  capacitor  modulation  and  through  inductor 


- 1.  ~  i - /n_r  n  o*  /  a  »r*i - i - 1 - - - - j..i  ...I  ... - -  ..i  *  t_ 

luuuuiaLLun  v rvtr-L  l.  «  ->  /  /  /  «  .Lilt:  luuui;lui  iuuuuici c i_uil  tJ-Ml L  U- 

its  an  almost  negligible  transient  if  the  frequency  shift  is 
small  in  comparison  to  the  resonant  frequency.  The  system 
whose  spectrum  Is  sought  uses  induction  modulation!  however, 
it  can  achieve  frequency  shifts  which  are  significant  in 
comparison  to  the  resonant  frequency. 


A  General  Switched-L  Circuit  for  FSK  Transmission 

The  circuit  of  Fig,  4  is  capable  of  a  zero-transient 
frequency  shift  for  the  cos(wt)  form,  or  those  denoted  by 
**B"  on  Fig.  3.  For  the  sin(wt)-type  frequency  shift,  or 
those  which  are  denoted  by  "A",  a  transient  response  will 
occur.  For  the  type-B  (cos(wt))  shifts,  the  initial 
condition  generator  must  create  a  current,  IQ,  in  the 
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switching  inductor,  L,  and  the  voltage  sources  must 

be  changed  simultaneously  with  the  L  switches.  1‘he  ad- 

s 

dition  of  L  will  lower  the  resonant  frequency.  Xo  return 
to  the  original  frequency,  the  ganged  switches  must  be 
placed  to  their  original  position  when  the  current  is  at  a 
positive  maximum.  However,  the  type-A  (sin(wt))  shifts  can 
be  achieved  if  the  current  source  is  a  short  circuit,  and 
if  the  frequency  transition  occurs  at  an  antenna-current 
zero-crossing  with  positive  slope.  Ihe  addition  of  L  will 

9 

again  lower  the  resonant  frequency,  and  the  transient  is 
caused  by  the  mismatch  in  the  capacitor  voltage.  The 
capacitor  is  charged  to  a  maximum  energy  level  at  one  fre¬ 
quency?  however,  the  frequency  change  corresponds  to  a 
different  capacitive  reactance  and  therefore  to  different 
voltage  and  energy  levels  from  those  which  are  necessary  for 
a  zero -transient  response. 
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In  both  the  example  cases,  the  voltage  sources  are 
assumed  to  have  phase  continuity.  That  is,  at  the  time  of 
each  frequency  change,  there  is  no  instantaneous  jump  In 
the  voltage  phase.  Such  a  constraint  implies  a  special 
synchronization  between  the  sources,  as  well  as  between  the 
circuit  current  and  the  switching  time  of  L,.  Appropriate 
synchronization  and  control  circuits  were  developed  for  the 
ETAM  (Ref  4i7). 

An  Antenna  Modulator  for  Synchronous  Rer.onation 

A  model  of  a  synchronously  resonated  antenna  modulator 
for  FSK  transmission  is  shown  in  rig,  5.  The  model  is  a 
variation  of  the  circuit,  in  Fig.  4,  in  that  the  effective 
primary  Inductance  is  if  the  switch  is  open,  and  the  in- 
ductsnce  changes  to  approximately  L ^  -  (M  /L2),  with  R2  very 
small,  when  the  switch  is  closed.  Not  shown  in  Fig.  5  is  a 


Rj  total  transmitter  and  antenna  resistances 

R2  resistance  of  Che  turn ng  coil 

C  antenna  static  capacitance 

L_  antenna  apparent  inductance 

Lj  sum  of  the  helix  and  tuner  inductances 

M  mutual  coupling  inductance 

L2  tuning  coil  self -inductance 

Fig.  5,  An  Induccively-Ccupled  Synchronous  Resonator 
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timing  circuit:  which  enables  the  switch  to  open  or  to  cluse 
and  which  changes  the  frequency  of  the  source  at  the  Instant 
of  an  antenna -current  zero-crossing. 

Operational,  synchronously  resonated  Ei>K  systems  u»c 
a  saturable  reactor  or  an  electronl<  ally  switched  inductive 
shorting  ring  technique  to  achieve  broadbunding.  Either 
technique  can  be  modeled  as  the  circuit  of  Fig.  5.  Such 
broadbanding  Methods  are.  however,  limited  by  the  large 
Inherent  saturable-reactor  time-constant.,  or  they  are 
limited  by  the  speed  of  the  solid  state  switches.  The  two 
time  delays  cause  an  error  in  the  switching  Instant  and  a 
prominent  phase  transient  occurs  for  transmission  rates 
higher  than  50  baud  (Ref  3i555). 

The  ETAM  has  demonstrated  that  VI. F  antennas  can  be 


synchronously  resonated  for  broad bunding  by  Inductors  which 
are  magnetically  coupled  to  the  antenna  tuning  inductance, 
and  that  are  short-  or  open-circuited  with  newly-developed, 
high-speed,  silicon-controlled  rectifiers.  The  result  is 
that  the  resonant  frequency  of  the  system  instantaneously 
follows  the  frequency  to  be  transmitted,  although  a  tran¬ 
sient  does  occur  due  to  the  capacitor  voltage  error.  The 
high-current,  high-voltage  rectifiers  are  triggered  so  that 
in  this  switched-L  conf iguration  all  changes  in  resonant 
frequency  occur  at  the  instant  of  an  antenna -current  zero- 
crossing. 
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Inipl  icntiona 

The  implications  of  electronic  broadband i  ng  are  that 
VLF  FSK  transmissions  can  now  bo  made  at  rates  much  higher 
than  previously  feasible,  and  in  addition,  the  shitt  between 
the  marking  and  spacing  frequencies  can  be  many  times  the 
system  bandwidth.  The  energy  spectrum  for  the  synchronously 
resonated,  FSK,  antenna  current  at  the  higher  baud  rate  will 
therefore  be  far  wider  than  the  bandwidth  for  an  antenna 
without  synchronous  resonation. 
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III,  Current  Analysis  for  Synchronously  Resonated  Antennas 


The  preceding  chapter  contains  a  circuit  diagram  for  an 
inductively-coup*ed  antenna  resonator  for  electronically- 
broadbanded  FSK  transmission.  The  purpose  of  Chapter  III  is 
to  develop  an  expression  for  the  antenna -current  transient 
response  of  antennas  which  use  the  resonator,  and  to  develop 
an  expression  for  the  capacitor  voltage  error. 

Antenna  Transient  Analysis 

If  the  antenna  circuit  is  considered  as  an  RLC  circuit, 
then  its  differential  equation  is  that  of  Eqs  (4e.)  and  (4b) . 
Gamble  (Ref  3iAppendi.x  B)  has  shown  that  if  the  steady-state 
current  is  v(t)  =  A  sin(wot),  then  the  steady-state  capac¬ 
itor  charge  is  expressed  by 

-  cos(v,0c)  (11> 
if 

<lr>2 «  ic  (12> 

Fig,  4  indicates  that  an  inductance,  L  ,  is  added  to  or 
removed  from  the  circuit  at  each  frequency  cha,  ge.  If  Eq  (4b) 
is  expressed  as 

q  +  2aq  +  w^q  »  ~  v(t)  (13) 

where 

o  =  ^  (14a) 
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( 14b) 


and  if  an  inductance,  L  ,  is  added  to  the  circuit,  then 
Eqs  (14a)  and  (14b)  are  changed  to 


al  = 


wi  =  Tl+l^c  (1 

If  the  source  is  simultaneously  set  to  v(t)  -  A  sin(wjt), 
then  Eq  (13)  is  changed  to 

q  +  2otjq  +  w^q  «  sin(wjt)  (] 

5 

and  the  Laplace  transform  of  Eq  (17)  is 

ot  2 

s  Q(s)  -  sq(o)  -  q(o)  +  2ajsQ(s)  -  2a^q(o)  +  w|Q(s)  = 

A;  7%  «' 

The  proper  capacitor  charge  for  a  zero-transient  response 
at  radian-frequency  is 


J |  sin(wjt)  dt  «=  q(o) 


in  accordance  with  Eq  (9),  However,  a  voltage  error  does 
exist  and  departure  from  Gamble’s  analysis  (Ref  3jAppendix  B) 
allows  the  initial  charge  to  be  redefined  as 


<l(o)  -  jg-U+Yj) 
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where  Is  the  capacitor-chargc  error-fraction.  The 
subscript  1  is  for  parameters  at  the  shift  to  any  general 
frequency  w^.  If  Eq  (15)  is  rewritten 

2  A  R  _  2A  _  i 

R  Tl+lJ  7  TT  al  K‘ 

5 

then  substitution  of  Eqs  (20)  and  (21)  into  (18)  yields 
Q(s)[s2  +  2a1s  +  w2]  r  g-j-d  +  YpCs  +  2aj]  - 


s  +  wi 


where  q(o)  **  i(o)  =  0  for  an  antenna -current  aero-crossing. 


Then 


Q(s) 


(s  +  aj)  +  aj 


*1  sZ  +  w?  1  (s  +  ct.)2  -  w 


1  (23) 


where 


w2  «  w?  -  a2 


The  inverse  Laplace  transform  yields  q(t),  and  differen¬ 
tiation  of  q(t)  gives  the  current  i(t).  Since,  from  Eq  (12), 

a*  «  Vj  (25) 


then 


wn  «  wL 


The  complete  response,  i(t),  is 


i(t)  »  [1  *  exp(-a1t)]  sintwjt) 

X 

where  Eqs  (15),  (16)  and  (20)  define  the  parameters. 
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'ii 


A  more  meaningful  form  of  Eq  (27)  is  produced  when  the 
relation 

B*  =  Ban<lwldth(3dB)  C28) 

<tt> 

is  substituted  into  the  exponential  term  of  Eq  (27),  The 
result  is 


i(t)  c  Cl  +  Yi  exp(-  tt  BW  t)]  sin(w^t) 


(29) 


Eq  (29)  emphasizes  that  the  transient  term  decays  more 
rapidly  when  the  antenna  bandwidth  is  increased. 

The  antenna  current  expressions,  Eqs  (27)  and  (29),  are 
based  upon  having  a  knowledge  of  the  error  fraction,  y^»  at 
the  time  of  the  frequency  shift.  Given  such  information,  an 
error  fraction  can  be  calculated  for  a  frequency  shift  which 
occurs  at  any  ensuing  antenna-current,  zero-crossing. 


Error  Fraction  Calculations 


The  capacitor  voltage  at  time  t  can  be  expressed  as 

vc(t)  =  g  J Q  i(a)  da  +  k  (30) 

where  k  is  the  voltage  at  t  e  0  and  l(t)  is  from  Eq  (27), 

The  constraints  cf  Eqs  (12)  and  (25),  with  the  requirement 

that  all  error-fraction  calculations  need  be  made  only  at  the 

instant  of  an  antenna -current  zero-crossing  (that  is,  when 

an  integer-multiple  of  cycles  has  been  completed),  yield 

*  2irna, 

**  Y_  a-ti  Cl  -  exp(- 

2?m 


t*=r 


a  CRawa 


-)]  ♦  k 


(31) 
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Fig.  6  depicts  the  time  relations  for  Eqs  (32)  and  (33) 
below.  From  Eq  (20),  the  initial  capacitor  voltage  at 
time  t  =  0  is 


waRaC 


*  (1  +  Y«) 


and  the  voltage  at  the  instant  of  a  later  frequency  change 
may  be  similarly  expressed  as 

-  dh<1  +  V  (3 


General  Error-Fraction  Time  Relations 
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Substitution  of  Eqs  (32)  and  (33)  into  Eq  (31)  yields 

R.w.  2n-rTa_ 

"  IT  e*P<— w— >3  -  1  <w) 

,  a  a  a 

from  which  the  error  fraction  at  the  next  frequency  3hift 

can  tye  calculated.  Eqs  (27)  and  (34)  provide  a  method  of 

calculating  the  antenna  current,  given  a  specific  KSK  bit 

stream. 


/ 

/ 

/ 


( 

/ 


1 


i 


! 

i 

! 
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IV,  Spectral  Analysis  for  Synchronously  Resonated  Antennas 


In  Chapter  III,  expressions  for  the  antenna  current  and 
for  the  capacitor  voltage  error  were  developed.  In  Chapter 
IV  a  method  Is  developed  for  finding  the  energy  spectrum  of 
FSK.  signals  which  are  transmitted  through  a  synchronously 
resonated,  VLF  antenna  that  Is  retuned  to  the  marking  and 
spacing  frequency  at  the  instant  of  the  antenna -current 
zero -crossing. 

Examination  of  the  Signal 

The  signal  set  which  determines  what  is  to  be  trans¬ 
mitted  Is  a  bit  stream  of  marks  and  spaces,  or  ones  and 
zeros.  Such  a  binary  bit  stream  is  similar  to  the  output 
from  a  teletype  or  from  a  binary  encryption  device.  For 
this  study  the  bit  stream  is  assumed  to  be  a  series  of 
independent,  equi probable  ones  and  zeros. 

By  assuming  that  n  (the  number  of  cycles)  is  infinite 
for  Eq  (34),  one  can  calculate  the  maximum  error  fraction  at 
the  instant  of  a  frequency  change  which  occurs  after  steady - 
state  oscillation  has  been  reached.  Consider  that  for  a 
change  to  the  higher  frequency  Eq  (34)  gives  the  maximum 
error-fraction  as  +0.1,  and  for  a  change  to  the  lower  fre¬ 
quency  Eq  (34)  gives  the  maximum  error-fraction  as  -0,1, 
However,  for  the  transmission  of  a  random  bit-stream  the 
error-fraction  would  not  reach  those  maxima.  The  values 
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that  the  error -fractions  assume  are  distributed  in  an 
unspecified  manner  between  the  two  maxima.  For  example,  if 
an  antenna  with  a  3  dB  static  bandwidth  of  65  Hz  is  synchro¬ 
nously  resonated,  then  the  time  constant  shown  in  Eq  (29)  for 
the  decay  of  the  transient  would  be  (l/Cn&K))  «  0,005 
seconds.  The  transient  would  exist  for  approximately  0,025 
seconds,  and  it  would  fully  decay  only  for  transmission  rates 
less  than  A0  baud.  The  transient  essentially  never  fully 
decays  for  htgh-baud  systems  such  as  are  within  the  capa¬ 
bility  of  the  ETAM. 

Comparison  of  Eqs  (27)  and  (34)  with  Fig.  5  reveals  the 

following  constraints  upon  the  system.  At  the  time  of  change 

to  a  higher  frequency,  the  switch  SW  is  closed,  thereby 

reducing  the  total  inductance  in  the  system.  However,  the 

switch  and  coll  resistance,  reflected  bach  into  the 

primary  in  accordance  with  the  mutual  coupling  between  the 

coils.  The  result  is  that  the  inductance  is  decreased  by 

L  and  the  resistance  is  increased  by  the  reflected  value  of 
s 

1*2 ,  so  that  the  system  bandwidth  is  increased.  The  converse 
happens  when  the  lower  frequency  is  switched,  i.e,,  the 
system  bandwidth  is  decreased.  Examination  of  Eq  (27)  shows 
that  the  steady-state  current  for  the  higher  frequency  is 
smaller  than  that  for  the  lower  frequency,  because  the  total 
circuit  resistance  is  larger  for  the  higher  frequency  signal 
and  the  peak,  source  voltage,  A,  remains  constant.  As  an 
illustration,  assume  for  a  specific  transmitting  system  that 
the  total  resistance  is  increased  by  15%  for  the  higher 
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Antenna  Current  Envelope 


GE/EE/73A- 10 


The  data  shown  on  Fig.  7  do  not  indicate  that  each 
frequency  shift  occurs  only  at  the  instant  of  an  antenna- 
current  zero-crossing.  Fig.  8,  however,  illustrates  the 
current/time  relations  near  a  key-shift  point.  The  fre¬ 
quencies,  decay  rates,  and  steady-state  current  levels  are 
for  illustrative  purposes  only.  The  important  feature  of 
Fig,  8  is  that  the  frequency  shift  must  occur  within  one 
cycle  of  the  frequency  which  existed  prior  to  the  key  shift. 
The  change  could  occur  anywhere  from  the  instant  of  the  key 
shift  up  to  one  full  cycle  later,  depending  upon  the  phase 
of  the  antenna  current  at  the  instant  of  the  key  shift. 


Fig,  8,  Downshift  and  Upshift  Antenna  Current  Relations 
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The  Stochastic  Process 

If  the  expected  value  of  u  stochastic  process  Is 
constant  In  time,  and  If  the  autocorrelation  function  of  the 
process  Is  a  function  only  of  the  time  difference  (t^  -  t2), 
then  the  process  is  wide-sense  stationary  and  the  power  spec¬ 
trum  for  the  process  is  the  Fourier  transform  of  the  auto¬ 
correlation  function,  Bennett  and  Rice  (Ref  li2355-2385) 
used  those  concepts  to  develop  the  power  spectrum  for  a  dif¬ 
ferent  type  of  FSK  signal  from  that  which  has  been  described 
in  this  paper.  The  signal  which  Bennett  and  Rice  used  was 
phase -continuous  (as  is  the  EXAM  signal),  but  the  frequency 
changes  occurred  at  the  time  of  the  key-shiftj  i.e, ,  there 
was  no  delay  due  to  the  antenna-current  zero-crossing.  In 
addition,  their  signal  did  not  include  any  transient  anal¬ 
ysis,  but  their  simplifications  and  the  postulation  of  a 
uniform  density  for  the  signal’s  phase  at  a  frequency  shift 
resulted  in  an  analytic  expression  for  the  described  power 
spectrum.  One  important  result  of  Bennett  and  Rice's  work 
is  that  the  analytic  expires sion  showed  that  the  power  spec¬ 
trum  for  a  certain  class  of  FSK  signals  decreases  rapidly 
with  respect  to  frequencies  away  from  the  transmitting 
system's  center  frequency.  That  class  of  FSK  signals  has  a 
modulation  index  (M)  that  is  0.5,  where 


M 


Baud  Rate 


(35) 


and  where  fu  **  upper  transmitted  frequency  and  fj  *  lower 
transmitted  frequency. 
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If  Eq  (27)  is  expressed  in  a  form  to  show  one  segment 
of  one  realization  of  the  antenna  current's  stochastic 
process,  then 

i(t)  ■  [1  +  y  exp(-at)]  sln(wt) 

o  S  t  <  2an  (36) 

Where  y,  a,  n,  w(  and  R  are  realizations  of  random  variables. 
However,  probability  densities  may  be  assumed  only  for  R,  a, 
and  w  because  they  are  functions  of  the  underlying  random 
telegraph  signal  whose  density  function  is  known  (R,  a,  and 
w  assume  one  of  two  values).  The  random  varible  n  may  have 
a  binomial -type  density  function,  but  there  is  no  definite 
data  to  confirm  that.  Finally,  y's  density  function  can  only 
be  described  in  terms  of  its  range  of  valuess  i.e.,  the 
maxima  described  under  the  previous  cub-heading.  As  a  result, 
the  mean  and  autocorrelation  of  the  ETAM  process  cannot  be 
calculated  directly.  As  a  hypothesis,  however,  the  ETAM’s 
power  spectrum  could  be  very  similar  to  that  for  Bennett  and 
Rice's  work,  because  the  range  for  y  is  in  general  a  small 
fraction,  and  the  range  of  values  for  n  can  only  be  different 
by  a  fraction  of  a  cycle  when  compared  to  Bennett  and  Rice's 
FSK  process.  Since  that  hypothesis  is  only  an  estimate,  the 
author  chose  to  develop  a  general-case,  iterative  method 
for  describing  the  power  spectrum  for  ETAM-produced  signals 
for  any  general  value  of  modulation  index  M. 
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Solution  Method;  A  Finite -Time  Tower  Spectrum 

An  expression  relating  the  energy  of  a  signal  v(t) 
to  that  signal's  Fourier  transform  is 


/>>«.*  /:  |.«-.i 


dw 


r  -Hw 

/  [t» 

J  -00 


,e(f )  j  i-  (Ira(f ) K  ]  df  (37) 


where  E  is  defined  as  the  energy  dissipated  by  voltage  f(t) 
applied  across  a  one-ohra  resistor,  or  by  a  current  f(t) 
ttirough  a  one-ohm  resistor,  Eq  (37)  states  that  the  energy 
cf  a  signal  is  given  by  the  area  under  the  |f(w)|  ^  (the 
magnitude  squared)  curve,  integrated  with  respect  to  the  fre¬ 
quency  variable  f  »  w/2ir  (Ref  6il27),  The  power  spectrum 
can  be  interpreted  as  the  energy  per  unit  bandwidth  (in  Hz) 
contributed  by  frequency  components  at  frequency  f.  The 
units  of  the  energy  density  spectrum  are  Joules  per  Hz, 

The  Fourier  transform  for  a  finite  number  of  cycles  of 
the  antenna  current  (an  n-cycle  segment  of  one  realization  of 
the  random  process  i(t))  can  be  expressed  as 


F(w) 


(HP)uw1  +  HvlWl 

+  j  +  j  Hy1w1 


(D  -  A BP  *  ACG) 

+  c2 

(ACU  -  C  ■»  ABG) 

~7  +  7 


(38) 


where 

2r<na, 

A  *  exp(-  — -~)  (38a) 

W1 
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B 

C 

D 

E 

G 

H 


al + 

2  2 
w^  -  w 

(38b) 

2wa^ 

(38c) 

cos( 

2tuiw  \ 

W1 

(38d) 

w2  - 

W1 

v/2 

(38e) 

sin( 

A1 

R, 

2ttiiw  s 

W1 

(  38f ) 

(38g) 

and  where  w^  is  the  radian  frequency  tor  chat  specific  seg¬ 
ment.  All  the  other  subscripted  parameters  must  be  assigned 
appropriate  values  for  the  conditions  at  the  beginning  of  the 
segment.  The  definitions  of  each  of  the  values  are  as  pre¬ 
sented  for  those  of  Eq  (27). 

The  transform  of  Eq  (33)  can  be  adapted  to  generate  the 


•r%  nt  ^  am  a  a*-  r%  I"  i  M  'i  »-  A 

j.'uwci  urn  x  i-»x  u  jl  xaix 


*  I  rtA  ^  >\M«  A  A  A  4T  «-Ia  A  1*M  A 

l-  ”  C  vl  100  VA  WWO  KUVO4W1IU0 


Eq  (27)  by  application  of  the  time-shift  theorem.  The  time- 
shift  theorem  states  that  if  a  waveform  is  shifted  in  time, 
then  the  Fourier  transform  of  the  shifted  wave  can  be  found 
by  multiplying  the  known  transform  by  exp(-jwTd),  where  is 
the  time  to  which  the  waveform  is  shifted.  Hence,  the 
Fourier  transform  for  a  finite  number  of  sequential  waveforms 
like  Eg  (27)  can  be  found  by  adding  the  exp( - jwTd)F(w)  prod¬ 
uct  for  each  segment.  That  is  possible  because  the  Fourier 
transform  is  linear  and  superposition  holds. 

A  finite  series  of  marks  and  spaces  from  the  underlying 
random  process  gives  rise  to  a  finite  sequence  of  antenna  cur- 
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rent-  waveforms  t£e,  (27)}.  Therefore  Eq  (38),  along  with  the 
time-shift  theorem, can  be  used  to  calculate  the  Fourier 
transform ' s  real  and  imaginary  values  at  specific  frequencies. 
The  value  of  the  energy  density  spec t Kura  at  a  specific  fre¬ 
quency  can  ther.  be  calculated  by 

F(f.)  *  **  ££  Re,(f.)'/  +  (2  lm,(f.)r  (39) 

j*=l  J  j=l  J  1 

where  Rej(f^)  and  Iro-jCf^)  are  the  real  and  imaginary  values 
of  tha  Fourier  transform  of  trie  j-th  segment  at  frequency 
Eqs  (38)  and  (39)  are  discontinuous  at  the  frequency 
which  is  being  transmitted.  An  expression  (via  I. ’Hospital 's 
rule)  that  is  valid  at  the  discontinuities  is 


lira  F(w) 
w^-w^ 


*=  Fy^w. 


[a2  (1  -  A)]  _  jnnH 
a2(a2  +  4w2)  1 


-  jHy1w1 


f*  f  •%  a  S  T 

L-uti  -  fl'j 

a4- (a2  +  4w2) 


(40) 


where  the  parameters  are  defined  as  tor  Eq  (38). 

In  order  to  make  the  finite-time  energy  density  spec¬ 
trum  meaningful  for  a  comparative  analysis,  the  spectral 
components  must  be  normalized,  and  that  is  done  by  adjusting 
the  spectral  values  at  each  frequency  so  that  the  area  under 
the  energy  spectrum  is  one  joule,  then  integrated  across  the 
entire  frequency  domain.  One  possible  way  to  do  that  is  to 
divide  the  energy  density  spectrum  of  Eq  (39)  by  the  total 
energy.  The  energy  (E)  of  one  segment  £Eq  (27)}  of  the 
antenna  current  is 


K 

I 
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E  «  {—  -  -pi^a-D  -  (41) 

W1  0^  +  40.*  J  4a^  +  /<aiwi 


where  the  parameters  are  as  defined  for  Eq  (38),  The  total 
energy  is  the  sura  of  the  energies  for  each  segment  of  the 

series. 

Kola tlons  Between  the  Circu it  parameters 

The  method  outlined  in  the  preceding  paragraphs  pro¬ 
vides  a  way  to  generate  the  energy  density  spectrum  for  the 
antenna  current  resulting  from  a  finite-length  random  tele¬ 
graph  signal.  There  does  remain,  however,  a  set  of  unspec¬ 
ified  relations  between  the  parameters  of  practically  all  the 
preceding  equations  The  parameters  are  «-a  and  Rb»  the 
total  circuit  resistances  for  the  case  of  the  lower  trans- 

/  #•  \  _  .  j  *.i u  *  *•  oni  *  ri-pH  fr^nurricv 

mitted  trequency  wu  - -  - 

(f  ),  respectively.  The  3  db  bandwidth  and  the  total  induc¬ 
tance  for  the  lower  and  higher  transmitted  frequencies  have 
been  designated  BWa.  BWb,  Lfl,  and  Eb,  respectively.  The 
first  relationship  between  the  parameters  is 


fh  - 

r.  *=  r  i  +  ,oi5(~Vt'“)j 

~b  “*  brv„ 


(42) 


The  relation  for  Eq  (42)  is  entirely  empirical  in  that  it 
was  formulated  from  the  results  of  Hartley's  work  which 
included  the  resistance  of  the  SCR  switches  (Ref  4*557). 

In  addition  to  Eq  (42),  Eq  (28)  can  be  used  to  express  the 

ratio  of  La  ana  Lb  as 
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a  Ab 


VWa 

VS 
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BW, 


fb  Rb  Bwa 
Ea 


(43) 


(44) 


The  use  of  Eq  (42),  (43)  and  (44)  provides  a  rational  basis 
for  selection  of  the  parameters  for  Eqs  (27),  (34),  (39), 
(40)  and  (41), 


Generation  and  Compilation  of  Spectra 

A  computer  program  for  calculating  the  energy  density 
spectrum  generated  by  the  transmission  of  a  finite-length 
random  bit-stream  through  a  synchronously  resonated  VLF 
antenna  was  written  and  used  to  provide  data  from  which  the 
spectra  for  various  cases  were  plotted.  The  principle  ex¬ 
pressions  used  in  the  computer  program  were  Eq  (38),  the 
Fourier  transform!  Eq  (39),  the  power  spectrumj  Eq  (41),  the 
energy  in  one  signal  segment i  and  Eqs  (42)  -  (44),  the  inter¬ 
parameter  relations.  The  program  is  compTetely  general  in 
that  the  modulation  index,  bandwidth,  baud  rate  and  the 
marking  and  spacing  frequencies  are  variables.  Appendix  A 
contains  a  listing  of  the  fundamental  portions  of  the 
program. 

Since  the  program  is  completely  general,  a  specific 
set  of  parameters  representing  those  of  a  hypothetical  VLF 
transmitting  system  were  chosen  to  demonstrate  the  program's 
spectrum-generating  capabilities.  The  hypothetical  system 
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was  designed  to  operate  at  27000  Hz  center  frequency  with  a 
3  dB  bandwidth  of  65  Hz.  Given  those  constraints,  Watt's 
work  (Ref  8*125-140)  was  examined  to  determine  a  reasonable 
value  for  the  antenna's  static  capacitance  and  the  value  cf 
0.05  microfarad  was  chosen.  Eq  (3)  then  gave  0.28  ohms  as 
the  total  circuit  resistance  at  the  lower  frequency.  The 
remaining  inter- parameter  relations  were  determined  in  the 
program  by  use  of  Eqs  (28),  (42)  and  (44). 

Trial  program  runs  showed  that  the  energy  density  spec¬ 
trum  for  individual  cases  reached  a  definite  pattern  after 
as  few  as  200  random  bits  were  read  and  processed.  The  spe¬ 
cific  cases  that  were  run  for  this  thesis  all  had  an  iden¬ 
tical  random  bit  stream  which  was  1000  bits  long.  That 
choice  provided  for  a  comparison  between  specific  cases  in 
which  the  baud  rate  and  the  modulation  index  {Eq  (35)}  were 
varied. 

The  energy  density  spectra  were  plotted  by  using  a  sep¬ 
arate  computer  program,  but  the  resulting  plots  were  not  well 
suited  to  comparative  purposes  because  the  spectrum  line  was 
Jagged  and  greater  than  one  inch  wide.  As  a  result,  a 
smoothing  routine  was  developed  for  the  plotting  program, 
and  that  particular  set  of  Fortran  instructions  is  listed 
in  Appendix  A.  The  entire  plotting  routine  was  not  in¬ 
cluded  because,  while  the  Fortran  language  is  general,  the 
auxiliary  plotters  are  certainly  not.  All  of  the  plots 
shown  in  the  remainder  of  this  thesis  have  been  smoothed 
by  that  routine. 


The  program  was  used  to  examine  several  cases  for  the 
above- specified  antenna  using  an  ETA.M-like  synchronous 
resonation  system.  For  the  first  set  of  six  plots,  Figs. 
9-14,  the  baud  rate  was  varied  from  125  to  2000  in  order 
to  examine  how  the  main,  or  center,  lobe  of  the  energy 
density  spectrum  varied  with  the  baud  rate.  The  next  three 
plots,  Figs.  15  -  17,  are  for  cases  where  the  baud  rate 
remains  constant  at  540,  but  the  modulation  index  is  varied. 

Each  of  the  plots  mentioned  above  is  a  portion  of  a 
normalized  energy  density  spectrum!  i.e.,  the  entire  spec¬ 
trum  has  one  joule  under  the  curve.  The  upper  and  lower 
transmitting  frequencies  are  marked  on  the  plots,  and  the 
normalized  energy  under  the  curve  on  each  plot  has  been 
integrated  and  entered  in  the  header  at  the  top  of  the  plots 
Due  to  the  finite  length  and  the  random  nature  of  the  bit 
stream,  the  integrated  energy  varies  from  0.92  to  1.04 
joules i  hence,  the  number  shown  on  the  plots  must  be  inter¬ 
preted  liberally.  The  "steps/shift  =  1"  entry  indicates 
that  the  frequency  shift  for  that  plot  was  an  instantaneous 
shift!  i,e.,  the  entire  shift  occurred  in  one  step  at  the 
instant  of  the  antenna -current  zero-crossing. 


□LOWER  FREO.  26865  HZ.  STPTIC  BRNDWIDTH  =65  HZ.  HOO.  INDEX  0.50 

X  UPPER  FREO.  27135  HZ.  GRFPHEO  ENERGY  =  1.0430  J.  STEPS/SHIFT:  1 


Fig.  10,  Spectra  for  540  Baud 


0  LOWER  FREQ.  2580 0  HZ.  STRT1C  SRNDWIQTH  =  55  HZ.  HOO.  INDEX 

X  UPPER  FREQ.  27200  hZ.  GRAPHED  ENERGY  =  Q.S21C  J.  STEPS/SHIFT: 
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25400  25 800  26200  26600  27000  27400  27800  28200  28Q00 

FREQUENCY,  HZ. 


CJ LOWER  FREQ.  266C0  HZ.  STATIC  BANDWIDTH  =55  HZ.  MOO.  INDEX 

X  UPPER  f RED.  27400  HZ.  GRAPHED  ENERGY  =  0.9974  J.  STEPS/SHIFT: 
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25400  25800  2620Q  26600  27000  27400  27300  26200  28600 

FREQUENCY*  HZ. 


Fig,  14,  Spectra  for  2000  Baud,  M  e  1/2. 


OWER  FRED.  2573G  HZ.  STATIC  BANDWIDTH  =  55  HZ.  HOD.  INDEX  1.00 

PPER  FREQ.  27270  HZ.  GRAPHED  ENERGY  =  0.7091  J.  STEPS/SHIET:  1 


FREQUENCY,  HZ 


□  LOWER  FREQ. ‘ 26450  HZ.  STATIC  BANDWIDTH  =  65  HZ.  HOD.  INDEX  2. CD 

X  UPPER  FREQ.  2'?540  HZ.  GRAPHED  ENERGY  =  0.6548  J.  STEPS/SHIFT:  1 
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Spectra  for  540  Baud 
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V,  Frequency  Modulation  for  Spectrum  Shaping 


A  method  for  determining  the  energy  density  spectrum 
for  the  transmission  of  FSK.  data  through  a  synchronously 
resonated  VLF  antenna  was  shown  in  Chapter  IV.  The  method 
was  programmed  in  Fortran  and  the  results  of  that  program 
appear  on  computer-generated  plots  under  the  last  sub¬ 
heading  of  Chapter  IV.  In  Chapter  V,  an  investigation  is 
made  to  determine  the  effect  of  frequency  modulation  on 
reducing  the  width  of  the  energy  spectrum's  main  lobe  and  on 
reducing  the  spectrum's  magnitude  at  frequencies  far  removed 
from  the  system's  center-frequency. 

What  Type  of  Modulation? 

The  sponsor  for  this  thesis,  RADC,  requested  that  the 
effects  of  making  the  frequency  changes  between  t'a  and  ffc 
over  a  discrete  time  period  (rather  that  the  total  shift 
occurring  at  the  instant  of  the  antenna -current  zero¬ 
crossing)  be  investigated  as  a  means  of  reducing  out-of-band 
spectra.  The  frequencies  f„  and  f^  are  the  lower  and  higher 
transmitted  signals. 

As  mentioned  in  Chapter  II,  synchronously  resonated 
antennas  could  have  a  theoretical  zero -transient  frequency 
change  once  each  cycle,  at  the  .  xt  antenna-current  zero¬ 
crossing.  Changes  that  do  not  occur  at  the  zero-crossing  can 
cause  significant  current  transiexits  which  reduce  the  signal's 
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detectability  and  potentially  can  cause  damage  to  the  SCR 
switches.  Hence,  a  feasible  and  desirable  method  of  fre¬ 
quency  modulation  for  a  synchronously  resonated  antenna  Is 
to  approximate  the  desired  frequency/time  relationship  with 
incremental  frequency  changes,  each  occurring  at  an  antenna - 
current  zero-crossing  and  each  lasting  at  least  one  complete 
cycle. 

There  are  limitations  as  to  how  much  frequency  modu¬ 
lation  can  occur  for  one  specific  key-shift,  because  for 
the  higher  key-shift  rates  at  VLF,  there  are  a  small  number 
of  individual  cycles  for  each  individual  baud.  For  example, 
a  1600  baud  keying  rate  transmitted  on  a  system  with  a 
center  frequency  of  27,000  Hr  results  in  each  mark  or  spa.ce 
having  16  cr  17  cycles.  The  frequency  modulation  certainly 
cannot  last  beyond  one  baud,  and  the  correlation  and  detec¬ 
tion  of  the  transmitted  signal  will  be  degraded  in  an  unspe¬ 
cified  manner  as  the  time  period  for  the  shaping  is  increased. 
There  is  simply  no  time  for  elaborate  frequency  shaping 
schemes. 

For  this  thesis,  therefore,  an  approximation  to  linear 
frequency  modulation  is  the  only  case  investigated  for 
spectra  reduction.  Fig,  IS  shows  a  frequency /time  chart  for 
(1)  a  segment  of  an  FSK  signal,  (2)  the  same  signal  with 
linear  frequency  modulation,  and  (3)  the  discrete  step 
frequency  modulation  that  is  within  the  capabilities  of  a 
slightly  modified  £TAM  circuit. 
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Fig.  18.  Multi-Step  Switching  for  FSK.  Signals 


Modification  of  the  Program 

The  program  was  modified  so  that  the  number  of  steps, 
or  the  number  of  discrete  frequency  changes  for  a  single  key- 
shift,  could  be  specified  by  a  data  card.  The  result  is  that 
the  finite-time  spectrum  for  a  random  bit- stream  could  be 
calculated  for  an  arbitrary  amount  of  frequency  modulation, 
or  an  arbitrary  number  of  discrete  steps.  The  transmitter 
frequency  shift  and  the  system  resonant  frequency  shift 
was  the  same  for  each  discrete  step,  as  shown  on  Fig.  18. 

The  spectra  for  several  cases  is  plotted  on  Figs.  19 
through  29.  In  each  case,  the  center  frequency  is  27,000 
Hz,  and  the  baud  rate  is  1600.  The  spectra  were  calculated 
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for  the  number  of  steps  per  shift  being  equal  to  1 ,  5,  10, 
and  15,  Three  frequency  bands  were  plotted  for  each  of  the 
four  different  cases,  the  bands  being  0  -  32000  Hz  (Low- 
Band),  25400  -  28600  Hz  (Mid-Band),  and  32000  -  64000  Hz 
(High-Band).  The  same  1000  bit  data  stream  that  was  used 
on  the  single-step  cases  was  used  on  the  multi-step  cases. 
In  addition,  the  error  fraction  was  calculated  for  each 
single  cycle  in  the  multi-step  frequency  transitions. 


0  LOWER  FREQ.  25500  HZ.  STATIC  BRNDWIOTh  =65  HZ.  MOO.  1N0EX  0.50 

X  UPPER  FP.EO.  27400  HZ.  GRAPHED  ENERGY  =  0.8935  J.  STEPS/SHIFT:  ] 


0  LOWER  FREQ.  26600  HZ.  5TPT1C  BRNOWIOTH  =  65  HZ.  MOO.  INDEX 

X  UPPER  FREQ.  27400  HZ.  GRPPHEO  ENERGY  =  0.8i74  J.  STEPS/SHIFT: 


□  LOWER  FREE.  25600  HZ.  STATIC  BANDWIDTH  -  55  HZ.  MOD.  INDEX  0.50 

X  UPPER  FREO.  27400  HZ.  GRAPHED  ENERGY  =  0.8205  J.  5TEP5/SHIFT:  ID 
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Fig.  23,  Mid-Band  Spectra  for  Steps 


□  LOWER  FREQ.  26500  HZ.  STATIC  6ANGW10TH  =  65  HZ.  MOO.  INDEX  0.50 

X  UPPER  FREQ.  27400  HZ.  GRAPHEO  ENERGY  =  1.0513  J.  STEPS/SHIFT:  10 
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□  LOWER  FREQ.  26600  HZ.  STATIC  BANDWIDTH  =  65  HZ.  HOD.  INDEX 

X  UPPER  FREQ.  27400  HZ.  GRAPHED  ENERGY  =  1.0262  J.  STEPS/SHIFT 
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0  LOWER  FREQ.  25600  HZ.  STATIC  BRNOWIOTH  -  65  HZ.  MOO.  INDEX  0.50 


FREQUENCY  *  HZ. 


0  LOWER  FREQ.  26600  HZ.  STATIC  BANDWIDTH  =  SS  HZ.  MOO.  INOEX 

X  UPPER  FREQ.  27400  HZ.  GRAPHED  ENERGY  =  0.0000  J.  STEPS/SHIFT: 


Fig.  27,  High-Sand  Spectra  for  Steps 


□  LOWER  FREQ.  26500  HZ.  STATIC  BANDWIDTH  =  65  HZ.  MOO.  INDEX  0.50 

X  UPPER  FREQ.  27100  HZ.  GRAPHED  ENERGY  =  0.0000  J.  5TEP5/5H3  FT :  10 
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VI.  Comparisons 


The  comparisons  are  divided  into  two  sections  which 
are  the  single-step  cases  and  the  multi-step  cases. 

The  Single -Step  Cases 

Figs.  9  through  14  are  plots  of  the  spectra  for  FSK 
signals  in  which  the  baud  rate  is  varied  and  the  modulation 
index  is  held  constant.  Figs.  15  through  17  are  plots  of 
FSK  spectra  wherein  the  modulation  index  is  varied  and  the 
baud  rate  is  held  constant.  Comparison  of  the  above  listed 
plots  with  the  equations  and  plots  contained  in  Bennett  and 
Rice's  paper  (Ref  l»2355-2385)  reveals  that  their  analytical 
expressions  (even  though  for  a  slightly  different  FSK  antenna 
current  than  that  for  the  ETAFl  synchronous  resonatlon  circuit) 
are  quite  adequate  to  express  the  general  characteristics  of 
the  spectra  lor  the  ETAM  circuit  for  the  cases  where  the 
frequency  shift  occurs  in  one  step.  Specifically,  the  main 
lobe  width  (for  the  cases  where  the  modulation  index  is  0.5) 
is  equal  to  1.5  times  the  baud  rate.  For  example,  Fig.  14 
is  a  spectrum  plot  for  a  baud  rate  of  1600,  arid  the  center 
lobe  is  between  25800  Hz  and  28200  Hz,  The  main  lobe  width 
therefore  is  2400  Hz  «=  (1600)  X  (1,5)  Hz.  In  addition,  the 
side  lobes  each  have  a  width  equal  to  0.5  times  the  baud 
ratet  i.e,,  the  sidclobes  for  the  spectrum  of  Fig.  12  are 
400  Hz  wide  while  the  baud  rate  is  800.  For  comparative 
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purposes,  a  plot  of  Bennett  and  Rice’s  equation  for  FSK 
spectra  with  a  modulation  index  =  0,5  is  shown  in  Appendix  B, 

Figs,  15  through  17  compare  closely  with  Bennett  and 
Rice's  figures  for  modulation  indices  of  1,0,  2.0,  and  3,0. 
The  impulses  on  the  spectra  were  very  narrow  in  the 
unsmoothed  data,  and  the  smoothing  routine  caused  a  slight 
offset  of  the  frequency  marker  on  the  figures  in  this  thesis. 

The  Multi-Step  Cases 

In  order  to  evaluate  the  plots  of  Figs.  19  through  29, 
an  additional  series  of  plots  was  prepared  on  which  the 
differences  of  the  spectra  for  two  plots  were  compared  in 
decibels.  The  spectrum  for  the  dividend  was  the  case  for 
Steps  *=  1,  and  the  divisor  spectrum  was  one  of  the  cases  for 
Steps  *=  5,  10  or  15.  At  a  specific  frequency,  if  the  spectra 
for  the  divisor  was  less  than  the  spectra  for  the  dividend, 
then  the  quotient  expressed  in  dB  was  positive.  Therefore, 
on  Figs,  30  through  38,  if  the  spectra  ratio  in  dB  is 
positive,  then  the  value  of  the  energy  density  spectrum  for 
the  multi-step  case  is  lower.  For  example,  Fig.  30  {the 
low-band  (0  -  32000  Hz)  comparison  of  the  Steps  e  1  case  to 
the  Steps  «=  5  case}  shows  that  at  22000  Hz  the  Steps  «  5 
spectrum  is  approximately  12  dB  lover  than  the  Steps  «=  1 
spectrum. 

An  evaluation  of  Figs,  30  through  36  with  the  FSK 
spectrum  in  Appendix  B  reveals  four  general  trends  for 
multi-step  frequency  switching.  The  first  is  that  the  main 
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lobe  ts  essentially  unchanged,  except  for  the  extreme  edges 
where  the  stepped  cases  are  as  much  as  seven  dB  lower.  Sec¬ 
ond,  the  first  five  sidelobes  on  each  side  of  the  main  lobe 
are  reduced  between  three  and  thirteen  dB  for  the  Steps  «=  5 
case,  between  six  and  eighteen  dB  for  the  Steps  »  10  case, 
and  between  nine  and  eighteen  dB  for  the  Steps  =  15  case. 
Third,  the  sidelobes  further  removed  from  the  center  fre¬ 
quency  than  the  first  five  sidelobes  are,  in  general,  reduced 
in  magnitude  by  multi-step  switching.  However,  not  every 
sidelobe  is  reduced  and  some  are  actually  increased.  The 
increases  are  indicated  by  spectra  ratios  that  are  negative. 
Fourth,  several  sidelobes  positioned  around  the  second 
harmonic  of  the  center  frequency  are  actually  increased  by 
ten  to  eighteen  dB  through  multi-step  switching. 


61 


... 


a  LOWER  FREQ.  25600  HZ.  gflUO  R.3TE-  1600  _ MOOULRTION  INOEX=  0.50 
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□  LOWER  EREGJ.  26600  HZ.  BAUD  RATE=  1600  MODULATION  INOEX=  0.50 


FREQUENCY.  HZ. 


lD  LOWER  FREQ.  25600  HZ.  BRUD  RSTE=  1600  MOOULRTION  1NDEX=  0.50 

X  UPPER  FREQ.  27400  HZ.  I  QUOTIENT  IN  OB  FOR  CSTEPS^  1  J  /  CSTEP5=  151 
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Fig.  38.  High-Sand  Comparison  for  Steps 
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VII .  Conclusions  and  Recommendations 


Electronic  broadbanding  by  means  of  synchronous  retuning 
is  the  only  practical  method  for  transmitting  high  data  rate, 
frequency  shift  key  signals  at  very  low  frequencies.  Effi¬ 
cient,  high  power,  high  data  rate  transmissions  are  not 
feasible  in  conventional  VLF  systems  because  of  a  severe 
transmitter/antenna  impedance  mismatch  and  because  of  pro¬ 
hibitive  signal  phase  distortion.  However,  because  the 
electronic  broadbanding  technique  shifts  the  antenna  system's 
resonant  frequency  in  accordance  with  the  transmitted  signal, 
the  bandwidth  of  the  transmitted  signal  will  be  wider  than 
the  system  bandwidth  for  an  antenna  that  is  not  synchronously 
retuned.-  This  thesis  contains  a  computer  program  method  for 
determining  the  energy  spectrum  foi  frequency  shift  keyed 
signals  transmitted  by  synchronously  resonated  VLF  antennas. 

Conclusions 

Single-Step  Cases,  The  approximation  for  the  energy 
density  spectrum  of  an  FSK.  signal  with  a  modulation  index  of 
0,5  that  is  transmitted  by  an  electronically  broadbanded 
antenna  has  shown  that  the  width  of  the  signal's  main  energy 
lobe  is  equal  to  1,5  times  the  baud  rate.  In  addition,  the 
sidelobes  each  have  a  width  equal  to  0,5  times  the  baud  rate, 
Those  spectrum  results  were  determined  by  an  iterative  method 
which  used  the  specific  stochastic  process  for  the  antenna 
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Electronic  broadbanding  by  means  of  synchronous  retuning 
is  the  only  practical  method  for  transmitting  high  data  rate, 
frequency  shift  key  signals  at  very  low  frequencies.  Effi¬ 
cient,  high  power,  high  data  rate  transmissions  are  not 
feasible  In  conventional  VLF  systems  because  of  a  severe 
transmitter /antenna  impedance  mismatch  and  because  of  pro¬ 
hibitive  signal  phase  distortion.  However,  because  the 
electronic  broadbauding  technique  shifts  the  antenna  system’s 
resonant  frequency  in  accordance  with  the  transmitted  signal, 
the  banawidth  of  the  transmitted  signal  will  bo  wider  than 
the  system  bandwidth  for  an  antenna  that  is  not  synchronously 
retuned.  This  thesis  contains  a  computer  program  method  for 
determining  the  energy  spectrum  for  frequency  shift,  keyed 
signals  transmitted  by  synchronously  resonated  VLF  antennas. 


Single-Step  Cases.  The  approximation  for  the  energy 
density  spectrum  of  an  FSK  signal  with  a  modulation  index  of 
0,5  that  is  t.?  isiritted  by  an  electronically  broadbanded 
antenna  has  shown  that  the  width  of  the  signal's  main  energy 
lebe  is  equal  to  1.5  times  the  bauo  rate.  In  addition,  the 
sidelobes  each  have  a  widtn  equal  to  0.5  times  the  baud  rate. 
Those  spectrum  results  were  determined  by  an  iterative  method 
which  used  the  specific  stochastic  process  for  the  antenna 


71 


GE/EE/73A-10 


current  of  a  synchronously  resonated  VLF  antenna,  in  which 
a  frequency  shift  could  occur  only  at  the  instant  of  an 
antenna -current  zero-crossing.  Also,  the  algorithm  included 
the  effects  of  the  current  transients  caused  by  an  inherent 
initial  condition  mismatch. 

Bennett  and  Rice  (Ref  li2384)  found  an  analytic  solution 
for  the  energy  spectrum  of  a  significantly  different  FSK. 
signal.  Their  analysis  did  not  include  any  time  delays 
caused  by  the  antenna-current  zero-crossing  criteria,  nor  did 
their  work  include  the  effects  of  the  antenna  current 
switching  transients.  However,  Bennett  and  Rice's  analytic 
solution  gives  the  same  width  for  the  main  and  side  lobes  as 
does  the  itera; ive  method,  and  the  values  of  the  two  spectra 
are  equal  within  one  dB.  Bennett  and  Rice's  analytic 
solution,  then,  provides  a  close  approximation  to  the 
iterative  method  used  for  this  thesis, 

Multi-Step  Cases.  All  irequency  shifts  for  the  single- 
step  cases  were  essentially  instantaneous.  The  multi-step 
cases  were  devised  to  determine  what  spectrum  changes  would 
occur  if  the  complete  frequency  shift  occupied  a  finite 
time  period.  The  program  was  modified  to  approximate  a 
linear  FM  frequency  shift  by  simulating  the  transmission  of 
a  series  of  single  cycles  at  incrementally  changing  fre¬ 
quencies  in  order  to  reach  a  different  FSK.  frequency.  The 
frequency  shift  increment  for  each  single  cycle  was  the 
3ame,  so  that  the  frequency  of  the  last  incremental  change 
was  the  terminal  frequency  for  that  shift.  For  example,  a 
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Steps  *=  3  shift  completed  one  cycle  at  the  original  fre¬ 
quency  plus  1/3  of  the  total  shift,  then  one  cycle  was  com¬ 
pleted  at  the  original  frequency  plus  2/3  of  the  total  shift, 
and  then  the  final  frequency  was  reached.  The  spectrum  for 
a  hypothetical  27000  Hz  center  frequency  antenna  which  used 
synchronous  retuning  to  transmit  FSK.  data  at  1600  baud  was 
calculated  for  Steps  =  5,  10,  and  15.  The  results  are  tnat 
the  main  lobe  is  essentially  unchanged,  but  the  first  five 
sidelobes  on  each  side  of  the  main  lobe  are  reduced  between 
three  and  thirteen  dB  for  the  Steps  =  5  case,  between  six 
and  eighteen  dB  for  the  Steps  10  case,  and  between  nine 
and  eighteen  dB  for  the  Steps  =  15  case.  The  other  side- 
lobes  are  generally  decreased  by  multi-step  switching, 
although  there  is  an  increase  in  the  spectra  for  all  the 
multi-step  cases  near  the  second  harmonic  of  the  center 
frequency  (54000  Hz),  The  first  sidelobes  on  either  side 
of  54000  Hz  were  increased  between  13  and  20  dB. 

Recommerdacions 

Examination  of  the  conclusions  and  criteria  for  this 
thesis  brings  forth  two  recommendations.  The  first  is  that 
the  detectability  of  the  signal  for  the  multi-step  cases 
should  be  investigated.  Even  though  there  are  significant 
Improvements  in  the  transmitted  spectrum  for  a  multi-step 
case,  the  approximated  linear  FM  in  that  signal  may  degrade 
its  detectability  by  standard  receivers.  Second,  the  por¬ 
tion  of  the  program  which  generates  the  linear  frequency 
shifts  could  be  modified  to  approximate  a  different  shift 
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pattern  so  that  further  improvements  in  the  energy  density 
spectrum  could  be  investigated.  Specifically,  the  sidelobe 
near  the  second  harmonic  could  possibly  be  reduced  by  a 
different  FM  pattern. 
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Appendix  A 

Computer  Programs 

The  main  computer  program  written  for  this  thesis  is 
listed  on  the  following  pages.  The  program  is  written  in 
Extended  FORTRAN  for  the  CDC  6600  computer,  hence,  some  of 
the  Input/Output  format  statements  may  not  work  on.  all  ma¬ 
chines  capable  of  compiling  a  FORTRAN  program.  The  type  of 
data  cards  required  are  defined  on  Comment  cards  within  the 
listing.  An  example  set  of  data  cards  is  shown  at  the  end 
of  this  Appendix. 

The  product  of  this  program  is  two  arrays,  one  for  a 
set  of  frequency  points  and  one  for  the  corresponding  set  of 
energy  density  spectrum  values  for  the  data  stream  fed  in  as 
data,  or  as  generated  from  a  random  number  function.  The 
author  did  not  include  an  output  format,  but  the  place  to 
insert  a  standard  textbook  variety  format  was  noted  in  the 
listing.  The  author  used  an  array  cnly  for  testing  of  the 
routine,  then  later,  he  used  a  Call  statement  to  a  subroutine 
for  creating  a  machine  plot  from  these  arrays. 

The  main  portion  of  the  smoothing  routine  for  the 
plotting  program  is  included  before  the  data  card  examples. 
The  spectrum  array  points  are  very  jagged  and  a  smoothing 
routine  is  necessary  for  a  machine  plot  or  for  a  set  of 
points  which  can  be  manually  plotted.  The  data  cards  shown 
for  the  program  contain  the  parameters  for  the  smoothing 
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routine,  hence,  the  parameters  are  available  for  the 
smoothing  routine  whether  it  is  used  in  the  main  program 
or  in  a  plotting  subroutine* 
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PPprPflH  POINTS  (INPUT, OUTPUT , TAPE  5- INPUT ,TAPE6  =  0UT°LT  , 

1HCT  , PUNCH, TAF'  ?,  TAPPET 

ARRAYS  3SA,  CL),  wCULcS,  AND  RATIC  MUST  PAVE  THE  SAME  LENGTH 
IN  THIS  °ROGPAP  AND  IN  T  H£  SUDNOL'T  INE  SPECTRA. 

COMMON  RcA(?0C3),  CLX ( 2  0  03)  ,  ENERGY 
C ] PENSION  JOULr S  (20 03)  ,  RATI  0  (  2 0  0  3) 
riMNSION  GAM  1(1),  CA(0(1>,  STPEAP(lCCn) 

LCL1VALENCE  (RLA(l),  JQULES(i)),  (CD(1),  RATIO(l)) 

INTEGER  3IT ,?ERCESf  CUES , OITO , STRE A" 

REAL  JOULES 

FIS  =  2, *3. 1415S2CS35 

INITIALIZE  INTERNAL  VALUES  BEFORE  PARAMETER  CARG  IS  REAC 
)  8  J  =  0 

FIT  TINE  =  0. 

TF  l  LAY  =  0 
l AST?  =  0 
LAST3  =  0 
7ERCES  =  0 

CNEE  =  0 

i  a  -  i* 

O  J  -  L 

K  s=  1 

ENERGY  =  0. 

REAC  200  ,  FH1GP  ,  F  l.  C  U  , A  MARK,  A  SPACE!  f6ALC»FNIN»FNAX,GAP'NAv 
1RN  ARK,RS°ACL,  A,ST£FS,PICT,F’^IN,SPCC1H,CEL,ERF00,RANCOP  , 
25NCTOP 

THE  PARAMETERS  ARE  OEFINEO  DElOH. 

FUCH  IS  THE  HIGHL*'  TRANSMITTED  FRFGLENCY 
F  l  (  K  JS  TH-Z  LCWeR  T  d  A  N  S  h  T  T  T  P  0  f-K’LGLfNCY 
ANAFK  IS  T.iL  ALPHA  PAT-V-tlTZP  F OP  T E  RICHER  FREQUENCY 
NHICH  IS  Afi-TIT  RARILY  CHOSEN  TC  Ft  TRANSMIT) EG  WHEN  A 
E INARY  ONE  IS  RlAU  IN  THE  PIT  STREAM, 

ASFACi  IT  THE  ALPHA  PARAMETER  FUR  THE  LGHER  FREQUENCY. 

FALC  ts  th~  PALO  R  A  T  t. . 

HMIN  IS  THE  LCNEri  F  CGE  FREQUENCY  CF  T^r  SPECTRUM  HThTCU. 

F  H  A  X  IS  THE.  UPPER  £  CGZ  FREQUENCY  CF  T  N r  CFRC1KUM  WINCCH. 

GAnHA  IS  the  INITIAL  ERRPP  pp  APT  I  Cl' » 

F  N  A  R  K  IS  THi.  T  C  T  A  L  P'STrTA'JCt'  FOR  T  F  u  HIGHER  FREQUENCY. 
FSPTCE  IS  THE  T  C  T  ?  L  R£S  I  STAN  J  £  f-  C  R  THE  ICNlO  FpEGU£NCY. 

A  IS  THE  TROnS'-ITTL  •  VOLTAGE  AmPlIUCE,  USUALLY  1.0 

SUFS  IS  THE  NL*'EC~  0"  FREQUENCY  CHANCES  1 0  COMPLETE  A 
SWITCH  PITH-’EN  FhIGH  AND  rLOW. 

nn  =  i  enables  a  plotting  su^pcliinl,  ttct  =  o  skips  it. 

r R 1 N  -  1  ENABLES  FRIhTING  of  PARAMF.T-KS  fcr  subroutine 

SPECTRA  CN  EACH  CALL  TC  IT.  9c  AN  ARE  CF  THE 
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VOLUME  OF  CAT  A  WHICH  MAY  PE  PRIMES.  IF  F  PIN  =  C, 

THFN  THE  PRINTEC  OUTPUT  IS  0MTT1EC. 

SFCCTH  =  i  fNAClES  SMOOTHING  IN  THE  PLOTTING  ROUTINE, 
f  F.  L  INVERSE  IS  THE  LENGTH  OF  THE  SUCTH1NC-  PAR. 

FRRCF  =  1  ENA  ELLS  c  R  FOR  FRACTION  CALCULATION  FOR  E  A  C  F- 
CALL  TO  SPECTRA,  ERROR  =  0  SETS  GAMMA  TO  ZERO. 

FAFCCM  =  1  ENAEl.ES  GENERATION  A  RANCC*'  E  IT  ST»EAW  PI 
CALLS  TO  TFE  RANDOM  NUMBER  GENERATOR.  RANCO^  =  0 
CAUSES  THE  r  T I  S  T c  E  A M  TO  "E  REAl  FRCM  CARCS. 

S»  CTOP  =  1  CAUSES  THc.  SMOOTHING  ROUTINE  TC  WCR*  ON  ALL 
THE  SPECTRAL  C.ATA,  SMCTOF  =  0  CALSES  THE  SHOO  THING 
ROUTINE  TO  SFCCTH  THE  VALUES  ANC  IGNORE  THE  PEAKS. 

**44*>**»-*‘****  #*«*  *♦**♦*#****# 

THE  ENO-OF-FILF  e L A G  FOR  TAPE(S)  IS  SET  TC  ONE  IF  THERE  IS 
NO  DATA  IN  THE  IN^UT  FILE. 

IF  »EOt:  (5)  «NE.  C  )  STOP 

PE  AC  201,  L,LS 

L  IS  THE  NUN0E='CF  FCTNTS  BETWEEN  FHIN  ANC  FFAX  AT 

WHICH  THE  S r EC TCUM  IS  CAUCULATEC.  L  MUST  NOT  L A R  C E c 
THAN  TH,-.  01HFNS1CN  OF  VHE  A  R  P  A  Y  PEA, 

LS  IS  THE  LENGTH  C  THE  FIT  ST RE  AN,  LS  MUSf  NCT  9£  GREATER 
THAN  THE  DIMENSION  OF  THE  ARRAY  STREAM, 

20  0  FCF> AT  (^1*5.0) 

201  F  C  c  F  A  T  (IS) 

C 

C  ESTABLISHES  R  A  "  I A  N  AND  FREQUENCY  F  A  P  A  ME  T  ER  SFT. 

WFARR  -  »I2*FUCH 
W  S  F  ACE  ~  ° I ? *  F LOW 
W  F  E  A  N  =  (WMAPK  +  HSFACl)  /?, 

RUN  =  (FMIN/WHE  AN)*PI2 

FFAX  =  (!~MAX/W  U  AN)  *P  12 

Wt'JF  =  RM  INV  KM  c  AN 

VMp>  =  RM  AX*  H  M  I  AN 

KELTA  =  (WHAX-WHIN  >  /  (u-i) 

FFfcflN  =  W  iCAK/f  1  ? 

FUN  =  WM IN/ PI  2 
FFAX  =  WMAX  /  R  I  2 
f CELIA  =  HOE  L  T  3  /  FI2 
C 

C  WMAFK  AND  AMARK  ARE  CALCULATED  INSTEAC  CF  FIXED  BY  WHAT 

C  WAS  READ  IN  PREVIOUSLY. 

FF  APK-R’S-’ACr*  (1 .  ♦  <F  HIGH- FLOW)  /  (2.  *  AS  FACS  /FT2)*.ni5) 
fF  AFK  = (PI  2/2, )  *  (  ( WMARR**2) »RMARK *  ( 2  *  A  SP  A  C  E / F  I  2 ) ) / ( ( W  S F  A  C l 
1 *« I)*RSPACE) 

C 

C  CL  C  ARS  ARRAYS  RF  A  AM)  Cl.X 

CC  bo?  MlUyL 
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F?/um  =  o. 

502  CLXHM)  =  0. 

PRIMS  TUT  THE  FARAmftERS  THAT  HlRL  READ  IN. 

FFIM  406 
406  FCRRAT  (Hi) 

PRIM  40  2,FhlGt- ,  ARARk',  RHARK,  FLOW  ,  ASFACE  »RSFACE» 
lECLC,FMAX,T/'MVA,A,L,FHIK,W'1IN1WMA)itFr:ELTA 
4C2  FORMAT  (H  ,T?0,*TH  r  A"3  AM  ETC  RS  R  E  AC  IN  HERr*  // 

IT  2  0  ,*  FHIGH*  ,  T  4  0  ,*  AM  ARK*  ,  T60,  *RrtARK*/ 
2Tic,3(1P*13.6,7X)  ,// 

3T l  C ,*FLOW* ,T4 C  ,  *  A  S  e  A<~E* , 760,*1SDACE*/ 

4T3c,3(iP£i3. 6 , 7  y )  ,// 

5T2C,*GAU1*,T4C,»FMAX*,T6  C ,  *  G  A  H*-’  A  *  / 
fcT 19,3(10-13.6 ,">X) // 

7  T  2  l,*AMOL  ITUnc  *  ,1  to  ,*  AROAY  LENGTH*, TE0,*FMIN*/ 
eTj°,iCF13.6,T2c,IE,TcT,  10E13.6// 
cTJp,*uMrj*,TO,*kMAy’',T6  0,*WDELTAV 
A T 1 S  ,3  (1P£  1  3.  6, ?>>////) 

TF  IS  FILLS  THJ  Ap  c  A  Y  STREAM 
622  IF  (R ANDO  1. EC. 0 . :  GC  T0  641 
LSM  =  LS  -  1 
CC  607  14=1, LS* 

F1P  =  IA 
F C  A  =  PA  (PI  A  ) 

FAN  =  RANF(POA) 

SIFFAM(IA)  =  1 

IF  (FAN. LE.  0.5)  STREAM(IA)  =  (i 
507  CONTINUE 

SIRIAM(LS)  =  3 
CC  TO  647 

641  FEAC  202,  (STREAM(IK) ,  IK  =1,LS> 

202  FCFRM  (7211) 

IF  (EOF ( 5  i  .NE. 0  5  Slur 

642  IF <FpiN.EQ. 0. >  GC  TO  607 
F  R  1  M  401 

409  FCRKAT(1U0,T6 , *STERS*,T 1 5,  *1 S*, T 2 3 , * HN*  , 

ITiE  ,  *  R  N*  ,  T  wr,  *  A  N  *  ,  T  5  r- ,  *  W  X*  , 

21f 5 ,*RX*, T7&. *  AX*  ,T56,*TCcLAY*, 

3TS5 ,*CNEN*,T 100  ,* ENERGY  *  ,/> 

PE  GIN  READING  AND  PROCESSING  THE  ST  t-_D  EIT  STPEAM 

TFt  FXAHINATICN  OF  TME  stre/m  ARRAY, 

607  FIT  =  STREAM ( K  > 

K  =  Ml 

IF  (MT. El.  1)  GC  TG  610 
IF  (EIT. El. 0)  IC  TC  61 3 
IF  (FIT, El. 2)  GC  TC  606 
IF  (EIT. El. 3)  0'  TO  609 

GC  TO  603 
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610  CNES  =  O'lES  *  1 
600  If  U.EQ. 1)  GO  TC  £03 
J  *  i 

FIIC  =  *IT 
GC  10  bO 7 


613  26  FCES  =  ZEROES  *  1 

GC  IP  6(10 

603  IF  (EIT.ET.9ITC>  GC  TO 
CC  IC  604 


607 


606 


C 

C 


IAST2  =  1 
GC  TO  604 


60S  l f  ST  3  =  1 


(^y)  VARIABLES  ARl.  FQp  SERIES  JUST  fi'.AC 
,fLn  VARIABLES  AFC  FOR  FPEVICUS  SEFf^S 
TC  EE  REA  0  FRC"  STREAM  NEXT, 


604  IF  (TITO.  ZQ.  1  )  GC  TO  611 
E1TC  =  GIT 

EIITI^L-  /wKUt'^b«ULJT  ii«iu 

7EFCES  -  0 
VNEV  =  mspace 

A  MV  =  AS  PALE 
FMV  =  R  S  °  A  C  E 
A C  L G  =  AM  A3K 
KLC  =  WMARK 
FCLC  =  RMARK 
GC  10  605 

C 

611  E  3  T  TTF£  =  ONES /E ACC* GIT T  IME 


GEES  = 
VMV  = 
AVEV  = 
FEtV  = 
PCLC  = 
FCLC  = 
FCU  = 
CC  TO 


0 

MM  A5  K 
A  M  A°  K 
RM  ADK 
AS ° ACE 
MS  PACE 
rspacz 

60G 


605 

614 


„'l  =  0  IF  THIS  IS  TH17  F I  PST  SERIES  TC  3 
IF  Ul.ro,  o>  GC  TC  c-14 
IF  (EPROR. Nl, 1 . 1  GNEW  s  0. 

CC  TO  623 
CMV  =  G  A  M  A 

]F  (ERROR. NZ« 1 «  >  GKt-F  =  0» 

J1  =  1 
GC  TO  023 
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IN  FROM  ST¬ 
OP  FCR  SERI E 


E  RE A C  F POM  ST R 


<n  in 


n  o  n  r»  r>  o 


Gh/lsb/  /D/Y-iU 


C 

623 


r 


C 

C 


r 


410 

C 

c 

421 


530 

C 

c 

633 

616 


C 

c 

c 

612 

C 

c 


c 

c 


CFANCH  OQINT  FCF  NULTI-STGP  SWITCHING. 

3  f  (STEFS.  “Q« 1 )  CC  TO  63’ 

istcp  =  srnrs  -  1, 

N  N  =  WOLD 
F  N  =  POLO 
n  -  AOLD 

CC  630  IS  =  1,  ISTEF 

CALCULATE  FRECL'ENCY  AND  RAMETERS  F  C  P  THIS  SINGLE  CYCLE  ANC 
ECR  NEXT  SINGLE  CYCL:-. 

TUS  LOO3  IN'OR  "  NENTS  THE  PARAMETERS  THROUGH 
(STEPS-1.)  CALLS  K  SPclCTRA  FOp  T  H  c  SINGLE  CYCLES. 

TFt  FARAMETFPS  PCR  tUIS  SINGLE  CYCLE. 

V  N  =  NN  +  U.  N  E  N -KCL  D) /ST  EPS 
AN  -  AN  +  (ANGW-ACLC) /STEPS 
FN  =  RN  +  (RNcW-FFLO) /ST EFS 

THE  PARAMETERS  FOR  TH£  NEXT  SINGLE  CYCLE. 


y> 

=  WN 

4- 

(WNPW-WCLC> /STEPS 

F  > 

=  RN 

4- 

(PNdW-PCLD) /STpPS 

A  > 

=  AN 

4- 

(  ANEW-ACLD) /STEPS 

CALL  SPECTRA  (  AN,WN,l,GNEW,L,TDELAY,V>MlN,WCtLTA,A,RN) 

TF  CFRTN.  EG.  0.  )  CC  TC  4?l 

FR1NT  413  ,STEF«  ,IS,  WM,=>H  ,AN,  WX,  RX,  AX  , TOE  LAY,  C  NEW  ,  F.NE  R  C- Y 
FCRNAT ( H  ,T3,lcE10.2,2X,I3,5X»lp93in.2) 

CALCULATES  GANHA  FOR  THE  NEXT  FPEClcNCY. 

CN  EW  =  ( (3X»WX  )  /  ( RN*  WN)  )  *  ( i.  +  G  NlH  -  cX  f  (  (- F  1 2  *  A  N  )  /WN)  )  -1 . 
IF  tERPOR, NE. 1 . )  CNEW  =  0 . 

TCELAY  =  TOE.LAY  +  FI2/HN 
CCNTINUE 

C6TERHINE  NUMEEC  CF  CYCLES  TO  MOVE  TCELAY  PAST  nITT IME  . 
F  =  1 

CYCLL  =  PI?/WN"W 

IF  (N^CYCLE  +  TOELAY  .GE.  ‘3ITTTME  )  C-0  TC  612 
N  =  N't-l 
CC  TO  616 


CALCULATION  OF  LCNG  SERIES  OF  CYCLtS. 

IF  (STlFS. EC. 1  ,  )  GC  TO  6  32 

CALCULATE  INC c  E  N  E  NT  AL  <N>  VALUES  ANC  KEEF  SANE  (X)  VALLES. 

K  N  =  WX 
AN  =  AX 
FN  =  PX 

THE  PARAMETERS  THAT  WILL  OE  HEEOEC  F  CR  THE  NEXT  SINGLE  CYCLE 
K>  =  WNFW- (WNEN-WCLC) /STEPS 
#>  =  ANEW- (ANSW-ACLC) /STEPS 
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C 

63? 


C 

63/* 


420 

635 


C 


R  >  ;  PNZH  -  (PNEV-FCLO/STPPS 
GC  10  634 

V>  =  WQLO 
F  )  =  R  OL  3 
A>  =  AQLO 
VF  =  VNEW 
FF  =  OFJEW 
A F  =  ANEW 
15  -  0 

CALL  SPECTRA  (AF,KN,N»GNrV,L,T02LAY,VMiFlfViCELTA,A,RN) 

IF  (FRIN.  SO.  0.  )  C-C  TP  635 

F  R  1  FT  410  ,CTEFS  ,1?,  WU,  RN  f  AN,  WX,  PX  ,A>  ,Tni_LAY,CNEW,  ENEPC  Y 
IF  (STEPS.  3(1. 1 .  )  GC  TO  o35 
FF1FT  42(1 
FCFFAT(1M  ) 

GFEV  =  ( ( RX*WX) / (RK*WN)  > Ml  .  +GMEW*EXF ( <-FI?*F*AN) /WF) )  -1 , 
IF  <FP°OR. NL. \ .)  CNEW  =  0  . 

TCclAY  =  TOLL  AY  +  N*nI2/WNEW 

IF  (LAST?.  SP.  1  )  GC  TO  619 
IF  (IAST3,  -a.l )  GC  TO  61<* 

GC  TO  o 0  7 


.CUES  CONTAINS  A  NOP  TAG  I7ED  ONE-STCEC  ENERGY  SFE  P  T  R  l  M 
T  F  A  T  IS  WORMALTZLC  10  ONE  JOULE  LFiER  T  P  L  CURVF.  OF 
JCILES  PER  l  Fill  FPcOUEMCY. 


fc IS  FRIFT  4tl,EN£PPY 

411  FC06T(// ,T5,  "TCT-U  STREAM  ENERGY  =  *  ,  1 F  E12 . 4 , 2X  ,  * 

1JCUF?.*//) 

CC  504  TJ  =  1  ,  L 

504  LIES (I J) =( <R - A(  I j ) **2+CLX< I J)  ^  ^ 2 ) /ENEFGY) *2. 


INTEGRATE  THl  NC°MALT2SD  ENERGY  SPECTRUM. 

AREA  =  0. 

IF  =  l-l 

CC  £06  IE  =  i,LM 

506  ARcA= (( C JOULES ( IE  +  11 +JOULESC IE) )  /2. ) *KOELTA) /FI2  +  AREA 


PRINT  412,APrA 

412  FORMAT  (14  ,Tr,*THE  ENERGY  IN  THE  NORMAL  I2EP  ENERGY  EAPC 
1  15  *,iP2i2. 4, 2X,  ’JOULES.*/) 

vCUFS  IS  SAVED  FOR  rOSSI0LE  PERMANENT  STORAGE. 

HIS  STORAGE  CIVICS  COULD  PE  A  MAGNETIC  TAPE  CR  AN 
CPERATING-SysTEM-poNTROLLED  FcRMfiNENT  CISC  FILE. 

FV  =  t.*AS',ACE/FI2 

VR  ITE  (  3)  L,FLCR  ,FhIGh,P  REAM,  STEPS  ,L  At  C,FM  IN',  F  MAX, 
15MCCTH,3ZL,T  drcc  ,LS,FW, smotop, ape  A 
VR  ITE  (3)  v  JLULFS  UT)  ,  IT  =  i,L> 

C 

IF  (RANOO 1.50. C.  )  GC  TO  643 
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C 

C 

C 


422 

643 

402 


413 


FINCHS  4  SIT  CF  CARDS  SC  THE  RANCCM  STREAM  CAN  RE 
IStO  AGAIN. 

FLNCH  (•??,  <STcEAM<IC>,  IC  =  1,LS) 

FCr|-AT  (7211) 

KC  =  K-l 

FRIM  40  3  ,  (SIR  :  AV(IZ)  ,TZ=l,KO) 

F  C  R  M  A  T  ( H  ,  T5,*the  aIT  STREAM  FOR  THIS  FCWF.R  SPECTRUM  *  , 
1J*NAS’//{T5,7211  )T 
FRIM  413 
FCKMAT  (1H-) 


60‘ 


IF  (FICT.ZO.O.  )  GC  TO  640 
THE  FILLS  T  N  E  ARRAY  RATIO 
CC  EOS  II.  =  1  ,L 

RtUC(IL)  =  (IL-1.)  *FOELTA  +  Fill  N 


WITH  FREQUENCY  VALUES. 


640 


FG3 


THE  IS  the  PLACE  an  a°fay  °r  int i k  g  routine  shculc 
EE  LOCATED  FOR  PRINTING  THE  CONTENTS  CF  THE  ARRAYS 
^CULE  3  AND  RATIO  (  TmL  SPECTpUM  ARC  THE  FRECUENCY 
VALUES)  THE  AUTHOR  ALWAYS  U^ER  A  FLCTIINC  °OUTIkE 
SO  AN  AP.PAY  Fl.CTTr R  IS  NOT  SHOWN,  THEY  AFE 
EASILY  ROUND, IN  FORTRAN  TEXT  B  C  C  K  S . 

YFIS  IS  THE  CALL  TC  THE  PLOTTING  RCLTINE . 

EMFC-Y  =  ARF  A 

CALL  FISTURE  <L  ,FLCW,FHIGh,FMEAN,  STEFS,SALC,FMN,FuAY,SKCTH, 
1CFL,  ERROR, LS,CM 

CLEARS  REA  AMC  CLX  ARRAYS 
CC  5  03  MN  =  1 ,L 
FFA(MN)  =  0, 

CLX'  <  M  N )  =  0. 


620 


FESET  INTERNAL 
EITTIHF  =  0- 
ENERGY  =  0. 

T  C  ELAY  *  0 
01  -  0 
0  =  0 

7ERCES  =  0 
CNCS  =  0 
M  =  1 
IF 


VALUES  DP10R  TO  PROCESSING  A  NnW  STREAM 


THE  LAST  BIT  IN  THE  STREAM  IS 
PROGRAM  LOCKS  FOR  ANOTHER  BIT 
THE  LAST  FI’  IN  THE  STREAM  IS 
FRCGRAM  locks  for  ANOTHER  SET 
A NOTH n R  St  T  CF  BIT  STREAM  CA°OS. 
THE  °E  ARF.  NC  CAFOS,  T  Etc.  FRCC-RAN 
(LAST2.  EQ.  1  )  CC  TO  620 
L  A  S  T  3  =  0 
CC  TO  60 3 
IASI?  =  0 
CC  TO  622 


IF 


IF 

IF 


A  TRC,  THEN  THE 
STREAM  CN  CARDS. 

A  THREE,  THEN  TH= 

CF  PARAMETER  CAROS  ANC 


STOPS, 
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v 


O 


c 


c 


o 

c 

c 

60  0 


c 

601 

c 


c 

c 

500 

c 

c 

c 


o 


m.EFOy  TI  N  ~  SFEr_T  ift  A NSW  ,  KNEW  ,  N  ,  G  C  L  t  ,  L  ,  T  C  1. 1  A  Y  , 

1 FF7V  ,  V.  P£  LT  A  ,  a  7  -b  i  M 
CCNNCN  REAtPOd?),  CLX(  2  0  03)f  energy 
REAL  IM 

FI?  =  2.*3.1AlFc2cc.3r, 

A 1  =  A 

F  =  Al/RNEw 

A  =  EXP( ( -PI ?*N*ANEK> /WNEW) 

rc  500  I  =  1,1 

V  =  KM  IN  +  (T-  1  )*KCLLTA 

P  =  ANEW**?  +  WN£K**2  -  H**2 

C  =  2.*ANEW*W 

C  =  (FI?*  N*W)  ,/WNFK 

f  =  KNEW*  *2  -  K**2 

F  =  COS( 1) 

G  =  ST  N  { 3 ) 

FIT  =  WNE  W- ( 1 . S-7) *WNEW 

LCLT  =  W  N  c  W  +  (l.E-7>*WNFW 

IF  (KjC-E.OUT.  ANE.W.LE.'JCUT)  GO  TO  SOC 

FE  =  (  (  (1.  -F)  *KNEK)/E*(GCLC*WNEW*  (?-A*E*F4A*C*G ))/(?**£ 
14C**?))*W 

IN  =  (  (  WNEW*G  )  /  •-+  (GOLP*WMLW*  (  A*C*F-C  +  f  »0*G)  )  /  (E**2+C**2)  )  *H 
CC  TO  601 


THE  IS  THE  LIHT  FUNCTION  WHEN  K  =  KNEW, 


h  \  =  Uj’a  n-'  »r  r*  •  ’  g  *  v  1  ♦  -  A; 


1  (A  NEW*  *2*  (ANlV**24E,*NN»K**2) ) > 

IN  =  <-PI?*N*H  /(?.  *KNEW)  ♦ 

1  (CCLO*WN"W*H*(-C*<l.-A)  )  )/(AfJCK**2* 

2  (ANEW  **?  +  4  .  *  KN‘:W**2)  ) 

GC  10  601 


Y  =  COS  ( W  *toe  l  A  Y  ) 
7  =  S1N(W*T0LL AY) 


FEA(l)  =  ( ~r*  Y  +  I M  *  ? )  +  REA ( I) 
CLMI)  =  (  I M * y  -  PE*Z)  +  CLX(I) 


CONTINUE 

CALCULATION  OF  THE  ENERGY  IN  THIS  SERIES  CF  NARKS  OF  SFACES 
C  *  (FI2*.N)  '  (  W N  •'  W *  2 «  ) 

F  =  (6,*G0LD*Kf  EN**2)/( ANEW**3  + 4 , *  A N E W * K N E K *  * 2) 

C  =  EXP<  (-PI2*N*ANEW)  /  (  W»  tLW)  )  -1, 

F  =  <  (GOLO*KNiK)  *  (GCL0*KNEH)  )  /  <  4  .  *  AN  t  K  *  *  3  . +«♦ ,  *ANE  W*  KN  E  K  *  *  2  . 
S  =  FXP( (-2, *F 1 2* ANEW*N ) /(KNEW) ) -1. 
p  N  E  FGY  =  ENEPGY  *  (h**2)  *(0-F*Q-F*S) 

A  =  A1 
FETIPN 
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S^OOTt'JNG  PC'JTrNr 

If  iSMcoTM.nr.  1 , )  re  TO  oM 

,RRAV  LI(IO)  ccktatns  Ix'wo^no  pcr. 

re  *23  M 3  =  1  » ? 

If  (f'P.rO,  1)  Gc  Tc  ry/tr 
I>  -  L/( 3 . "DLL  ) 

If(IX,LT,4)  G  0  TP  C5S 

If<lX»Gr«50)  GO  TC  63^ 

It  =  1 
It-  =  IX 
cc  Tc  fc>51 
‘  J>  =  3 
IL  =  1 
IS-  *  IX 

ICCC  ARR4Y  LT<lt), 

re  E2o  mx=i,ix 

*■  I  ( f-  X )  ~  R~A  (  1  L  -l+HX) 

LUJX)  =  LI  (IX) 

II<1)  =  LIC1) 

riCCLE  v  i  LUr^  ^  c  F  Llin^  VUu-s  FCR  LT(lC)  anc  test  Av;mm 

M?  =  lx- ?  • 

CC  5?J  MY  =  1,»; 
f'S  =  MY  4-  i 

7f  <SfiCTO^'-c!l.)  G0LT0I£2l"  LI(1,,)  7  (lx'n  *  MV 
ca^|,*GT*lT,^>)  r’°'TC  652 

CCCE22SMAMiDJt^?CLA  °  V£l-U'S  70  EC!-fl>-  MIDDLE  LTC)  VALLcS. 
MI  =  MA  4-  IL 
K  =  MA  4-  l 

FEA(MT)  =  LTC-C) 

CtMINuE 

Thr  FcmtRS  «**  t^t  FcK  fkc  er 

Ip  =  IH+1 

If  tlH,Li.  L)  GC  7C  6  5  0 
C  CM  I  MU* 

CCniKUZ 
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******  *  *  *  *  PARAMETER  CARDS  *********** 


PARAMETER  CARS'  F  CP  A  TYPICAL  PUS  A  R  c  S  F  C  F  N  FELOW. 
AMAKK  /NC  F  M  A  3  <  AP'  lL'FFY  ,'  A ,J  C  5  ONLY  IF  TFCSE  7  r  C  V  CLUE  S  A  F  E 
NOT  CAICILATCO  WITHIN  Tpr  P^G3Au,  NOTE  THAT  L  AN  F  LS  ARE 
INTEGER  VALUES.  TP;  c  I  T  STR.AP  CONTAIN*7  FEE  F  A  N  0  C  F  PITS 
WHICH  ACl  E.'NJEJ  -3  Y  A  (  )  TO  CAUSE  A  N  E  W  F  A  F  A  ;i  r  T  c  q  cry  ANO  A 
FEW  FIT  5  1  c  E  A  M  TO  RE  FEAR  IN.  IF  LESS  ThAN  LS  FITS  ARE  T  N 
YOUR  F  A  F  TITULAR  3IT  STREAM  ON  CARDS,  YCL  FUST  STILL  RcOVTCE 
FMCLCF  FLANK  CARTS  CR  PLANK  BITS  TO  FILL  LS  WCFCS  CF  THE 
AKPAY  STREAM. 


27  A  0  0  . 
2oC  C  3  . 
25  0. 

2  0  A .  1 
160  0  . 
25  A  0  0  . 


2060  F . 

0. 

2. 

.28 

1. 

1  . 

0. 

0. 

1  *  'v 
25 . 

1. 

C. 

1 . 

AO  1 
1000 


FFIC-T 

PICK 

AFAFK  DUMMY 

ACFACE 

E  ALC 

F  F  IN 

F*<AX 

G  A  F  ^  A 

F  v  A  F  K  DUNN  Y 
R  SPACE 
A  AFFLITIJOE 
STEPS 
RIFT 
F  c  1 N 

SvOCTh 

CEL 
£  E  RCR 
RANDOM 
S  M  0  T  C  F 
L  CAFE 

LS  STREAM  LFNJTH 


f 


column  .*ti 


column  **  7 Z 


Ir 


ii ni  c c 1 1  j ci  oil i  n n  1 1 1  inoa  cr 


iOinooi  ricijnoioii! 


oooi l  c  re  uoono mi  j  c Cirr 

moi  i  h  f  iccim  lullin' 
uiioi  n  i  ci  oi  u  uoo  i  ii' 

3101  r  C  0 1  CC1C011100000* 

oioocciicirotnoiooci 

11101 1 CC 1 C1P11 0 1  TOO  0  1 
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Appendix  13 

Ideal  Spectra  for  M  =  1/2 

The  plot  of  Fig.  B-l  Is  a  normalized  power  spectrum  for 
a  constant  amplitude  frequency  shift  key  signal  with  a  modu¬ 
lation  index  of  1/2.  The  plot  was  provided  by  John  Gamble, 
the  sponsor  for  this  thesis.  The  analytic  expression  for 
the  spectrum  is  referenced  as  Bennett  and  Rice’s  werk 
(Ref  Ii2384). 

The  plot  is  normalized  to  a  vaiue  of  one  watt  per  unit 
frequency  at  the  center  frequency,  F  .  That  normalization 
technique  is  different  than  that  used  for  the  plots  in  the 
body  of  this  thesis ;  they  were  normalized  to  one  joule  under 
the  entire  spectrum  of  joules/unit  frequency.  The  change 
from  joules  to  watts  makes  no  difference,  however,  because 
unit  power  is  the  integral  of  unit  energy  for  one  second; 
l.e, ,  the  normalized  spectrum  is  the  same  regardless  of 
whether  joules/unit  frequency  or  watts/unit  frequency  is 
platted  on  the  ordinate.  The  important  information  conveyed 
by  Fig.  B-l  is  that  the  main  power  density  or  energy  density 
lobe  width  is  1,5  times  the  baud  rate  (B) ,  and  that  the 
side  lobe  widths  are  0, 5  times  the  baud  rate. 
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